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THE  PROPAGATION  OF  ERRORS,  FLUCTUATIONS  AND  TOLERANCES 
BASIC  GENERALIZED  FORMULAS 

ABSTRACT 

^  The  classical  formulas  for  the  "propagation  of  error" 
permit  the  approximate  calculation  of  the  variance  of  a  function 
of  variables  whose  variability  is  known.  The  adequacy  of  its 
approximation  has  often  been  doubted.  Generalized  formulas  are 
here  obtained,  not  only  for  average  values  and  variances,  but 
also  for  the  next  two  curaulants. 

The  classical  formula  turns  out  to  be  better  than  expected 
--  and  further  approximations  can  be  confuted  when  necessary. 

Expressing  the  individual  variables  in  well  chosen  terms 
sinplifies  the  generalized  formulas  and  increases  the  accuracy  of 
the  classical  one. 

Higher  cocumulants  are  introduced  to  aid  the  algebraic 
work  of  derivation. 

k 

\ 


Abstract  -  2 


The  classical  formula  for  tne  "propagation  of  error", 
where  y  *  f(x1,  x2,  ...,  x^,  and  the  x&  suffer  Independent 
errors  or  fluctuations,  is 

var  y  £  2  (|:  -)  var  x  :  *  2  f  2  var  x 

where  the  subscript  o  means  i  hat  xg,  ...  xk  are  at  their 
average  values  and  where  we  ha  *e  written  f&  for  the  value  of  the 
a-th  partial  derivative  at  this  we  rage  point.  If  the  errors  or 
fluctuations  are  correlated/  this  generalizes  to 

var  y  -  Stf^  var  xa  +  2  oov  [x^j 

where  2*  indicates  summation  over  all  distinct  terms  with  ap¬ 
parently  different  subscripts  actually  different,  (see  Sec¬ 
tion  38  i or  further  details).  These  formulas  \re  exact  when  y 
is  a  linear  function  of  the  x‘s. 

7  "  ao  +  al  X1  +  x2  +  -•  +  ak  \ 

and  inexact  otherwise.  They  have  frequently  been  used  for  pre¬ 
dicting  the  variation  in  a  response  y  due  to  errors  or  to  wlthln- 
t  ole  ranee  fluctuations  in  the  x's.  There  has  been  a  tendency  to 
regard  their  use  with  suspicion  because  of  this  inexactness.  There 
has  been  no  indication  of  the  extent  of  inaccuracy  to  hold  this 
suspicion  in  proper  check. 


Attract  -  3 

C’he  present  memorandum  develops  formulas,  formally  cor¬ 
rect  through  tenns  of  order  rr> ,  for  the  average  value  o  1  y,  its 
variance,  ana  its  next  two  higher  cumul&nts.  The  most  important 
conclusion  is  that  the  classical  propagation  formula  is  much  better 
than  seems  to  be  usually  realized.  Examples  indicate  that  it  is 
quite  lixely  to  suffice  for  most  work.  The  generalised  formulas 
allow  us  (i)  to  check  the  accuracy  of  the  classical  formula,  and 
(ii)  to  obtain  much  more  accurate  results  in  the  few  cases  where 
they  may  be  needed. 

While  the  classical  fomula  is  often  quite  accurate 
when  used  in  the  terns  in  which  the  problem  originally  appears,  it 
is  often  possible  to  improve  its  accuracy  by  a  better  choice  of 
terms.  If  each  of  the  individual  or  component  variables  x.j ,  xg, 

...,  xk  is  expressed  in  terns  of  its  partial  effect  on  the  res¬ 
ponse,  or  in  such  terns  that  its  partial  effect  is  linear,  we' 
shall  say  that  the  terns  are  well-chosen.  Unless  the  design  value 
of  one  or  more  x's  has  been  chosen  to  make  y  a  maximum  or  a  mini¬ 
mum,  it  will  usually  be  possible  and  practical  to  use  well-chosen 
terms.  In  treating  the  delay  time  of  a  delay  line  made  up  of  LC- 
sectlons,  for  example,  we  are  led  to  express  component  variables 
in  terns  of  the  sqtiare-root  of  inductance  or  the  square -root  of 
capacitance.  In  treating  the  attenuation  of  a  ir-section  attenuator, 
for  further  example,  we  are  led  to  specify  the  shunt  elements  in 
terns  of  their  conductances  and  the  series  element  in  terns  of  its 
resistance. 


Ke’PI“«  furth,r  retire  .tension  to  the 

a kexmts  of  distribution  of  the  individual  ..„h,ku. 

vw  ±«<u.vxauai  variables,  measured 

conveniently  and  dlmensl onl ta sly  by 

ave  (xa  -  x  )3 

■V  «n»  .  £1  Q, 

a  i^i*; .  xa)^,^ 

Using  2*  again  to  indicate  sum-nation  over  all  distinct  tents 

with  apparently  different  subscripts  actually  different,  using 

7a  and  ca  to  describe  the  distribution  of  the  a-th  Individual 

variable,  and  abbreviating  the  various  partial  derivatives  of  f 

evaluated  -hen  ij .  ^  are  each  at  Its  average  value  as 

illustrated  by 

f>  _  (  d^f  ) 

aab  (3x^”3xa  ax^)  a11  x‘s  at  average  values, 

we  can  write  that  propagation  formula  for  the  variance  which 

applies  when  (1)  fluctuations  In  the  various  *a  are  Independent, 

and  (11)  the  *a  have  been  expressed  as  to  make  partial  effects 
linear  quite  simply  as  follows: 


var  y  *  2  f  2  a  2 

a  a 

+  ?*(faabfb  +  "ab"  +VanK2°b2 


Vaabb+  fab  faab  +  J  faaa bfb  >7^^ 


+  terms  of  order  >  a 


AfettHtct  •  5 


Various  aspect*  cf  this  fcmuia  Jestrve  attention  s 

(1;  rhe re  is  no  ter®  of  order  or.  Under  our  assump¬ 
tions  the  first  correction  ten.  is  two  order©  higher  than  the 
leading  tern.  Iftus  the  leading  terra  does  better  than  we  ©light 
expect . 


(2)  The  first  correction  ter®  involves  both  second 

(fab^  aaa  ^hlrd  (faab  and  fabb^  derivativ on  an  equal  footing 
and  with  equal  coefficients.  If  we  arc  to  go  beyond  first  deriva¬ 
tives,  we  should  go  to  the  third  derivatives. 

(3)  Through  terms  in  cr^,  which  is  usually  further 
than  we  have  any  excuse  to  go,  we  need  only  the  values  of  the 
a  and  the  r  .  Fewer  quantities  enter  than  we  might  fear. 

(4)  Much,  but  not  all,  of  this  simplicity  conies  from 
express  the  individual  variables  in  proper  terms.  The  general 
formula  would  have  8  terms  instead  of  3.  Correct  choice  of  terns 
is  very  helpful  but  not  essential. 

The  parallel  formula  for  the  average  is 
ave  y  =  yQ 

+  faabb>«aV 


+  S*(T2  faaabb)7-0' 


3.  2 
a  'a  ab 


+  terms  of  order  >  . 

Analogous  formulas  for  the  third  and  fourth  cumulants  are  given 
in  the  body  of  the  memorandum.  (Table  1,  page  II-6), 


«  t 

It  U  oft**  i-  'h  eliminate  the  fr.rtt 

correction  term  in  the  . '?•  t,zi  via  by  changing  tne  terns 
in  which  the  response  is  express?!.  (In  the  examples  cited  this 
leads  to  multiples  of  (1)  the  cu^e  of  a  constant  less  the  delay 
time  and  (li)  the  cube  root  of  tl  e  attenuation  plus  a  constant.) 
While  the  possible  advantages  deserve  being  kept  in  mind,  the 
first  correction  term  is  usually  small  enough  so  that  this  re¬ 
finement  is  unnecessary.  (The  general  question  of  simplifications 
is  discussed  in  Section  9,  pages  II.  9-11.) 

If  the  Individual  variables  are  not  necessarily  ex¬ 
pressed  so  as  to  make  f  ,  f  etc.  vanish,  but  do  have  inde- 
pendent  fluctuations,  the  propagation  formulas  are  less  simple, 
but  can  be  expressed  as  in  Table  A,  following.  (Similar  for¬ 
mulas  for  propagation  into  higher  cumulants  are  given  on  page 
1-8.)  These  formulas  require  more  information  about  the  fluctua¬ 
tions  of  the  individual  variables  —  information  provided  dlmen- 
sionlessly  by  the  values  of  F  and  Q  ,  where 

a  a 
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Propagation  formulas  for  the  case  where  the  fluctua¬ 
tions  of  the  individual  variables  are  not  independent  are  of  in¬ 
creased  complexity ,  and  will  be  reported  later* 

The  c  mputations  leading  to  the  formulas  are  simpli- 
lied  by  the  introduction  for  the  next  two  higher  cumulants  of 
fco-’  quantities  related  to  thr-e  cumulants  as  the  covariance 
is  related  to  the  variance*  '/hay  involve  3  and  4  arguments, 

respectively.  (Definitions  may  be  found  in  Section  19,  pages  V. 

1-2). 


The  structure  of  the  detailed  account  is  indicated  by 
the  following  table  of  contents.  Results  are  stated  and  dis¬ 
cussed  in  the  first  four  parts,  with  details  left  to  the  remain¬ 
ing  parts. 

The  examples  have  bear,  treated  both  by  "main  strength 
and  awkwardness",  and  by  more  polled  methods,  so  as  to  give  some 
idea  both  of  what  would  be  required  for  less  simple  examples  and 
of  how  easily  examples  can  be  handled. 
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1.  Introduction 


Fluctuation*  in  "individual”  or  "component"  variables  often 
combine  to  produce  fluctuations  In  a  ''resultant",  "system"  or  "overall" 
variable  --  fluctuation*  In  a  response.  Two  classical  examples  are 
the  combination  '/propagation")  of  errors  of  individual  operations 
in  a  physical  measurement  and  the  combination  of  within -tolerance 
fluctuations  in  a  mechai leal  or  electrical  assembly.  In  these  situa¬ 
tions,  as  in  a  wide  variety  of  others,  it  is  ofte ■  desired  to  relate 
the  distributions  of  component  fluctuations  to  the  distribution  of 
resultant  fluctuations.  For  a  long  time,  physics  has  used  the  formulas 
for  "the  propagation  of  error".  With  the  rise  of  modern  statistical 
theory,  it  has  been  natural  to  regain  the  propagation -of -error  formulas 
as  doubtful  first  approximations,  and  to  try  to  avoid  their  use  in 
situations  of  complex  dependence,  even  at  the  cost  of  proposing  great 
complications  in  stuty  or  experiment  or  even,  perhaps,  at  the  cost  of 
seeking  no  answer  to  the  problem. 

Insufficient  attention  seems  to  have  been  gi1  en  to  the  <jues- 
tion  of  how  accurate,  or  inaccurate,  the  propagation -of -error  formulas 
would  be,  or  tc  how  tney  may  be  improved.  The  obvious 


-  !  -I  * 


first  •tifp  is  to  ofetain  ar.i 
thin  is  acne,  the  classical 
most  of  us  Hal  suspect*  i. 


utttdy  better  approxiamtions.  men 
formula*  turn  out  to  b*  better  than 


But  more  can  be  done.  The  typical  problem  is  naturally 
expressed  as  follows; 


Given  the  functional  relation 


z 


*  *  *  > 


i 


and  given  Information  about  the  distributions  of  the  Individual 
w's  which  distributions  have  been  so  chosen  as  to  make  the  values 
of  z  be  distributed  near  zQ  what  Is  the  probability  that,  given 
a  "tolerance"  +  6  [or  -e  or  both],  that  z  will  exceed  z  +  B  [or 
fall  below  zQ  -  e,  or  either]? 

There  Is  nothing  in  the  expression  of  this  problem 
which  requires  us  to  work  in  terms  of  z,  rather  than  z2,  log  Zj 
v/T  or  z3:  We  can  translate  the  boundary  value  z  +  6  for  z  -  ej 
into  each  of  these  scales,  or  any  other.  Similarly,  we  may  re¬ 
place  W;L  by  w22  log  wx,  or  any  other  helpful  function.  And  we 
may  do  the  same  with  wg,  w y  . . „ ,  wR.  Thus  we  may  as  well  work 
with  variables  y,  uk,  a  functional  relation  y  =  f(u  , 

u2,  uk),  and  a  tolerance  yQ  +  6*  for  yQ  -  e*]  if  we  gain  by 

doing  this.  It  will  appear  that  we  can  increase  the  quality  of 

approximation  of  the  classical  propagation -of -error  formula  in 
this  way. 


■16 : 
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‘A'hile  the  same  techniques,  concepts  and  Insights  could 
be  app^iea  to  the  combination  of  (statistically)  dependent  fluc¬ 
tuations,  the  problem  of  combining  f statist  ?  caliy)  independent 
fluctuations,  --  which  leads  to  sir..,,j  er  results  --  is  of  suffi¬ 
cient  importance  to  lead  us  to  treat  only  that  case.  (After  all, 
most  practical  problems  with  depe^  .  nv  fluctuations  can  be  con¬ 
verted  into  problems  with  indecent  -  .  .,nes  by  pushing  back  to 
suitable  variables  --  often  to  varied  os  earlier  determined.) 

I  PROPACM^ZV.. 


2.  Moments 

The  basic  tool  of  "proper-  ..  of  error"  has  always 
been  the  "mean  square  deviation",  ci  .  as  statisticians  now  say, 
the  variance".  If  we  are  to  go  tn  ,  ^re  detailed  approximations, 
we  must  expect  to  supplement  this  wi r~i  more  complex  quantities. 

The  analogy  between  probait  _ty  distributions  and  mass 
distributions  and  the  moments  of  for.-,  area,  etc.  of  mechanics, 
early  led  to  the  description  of  prob-t.Uty  distributions  in 
terms  of  moments* ,  either  about  zero 


t 


i 

sp(z)  dz 

L 


'\  vre 


22  C  rp 

v  « 


z 


2 


fz^p(z)  dz  «  avj 


*  frnrrrfrrfiTiti^  nuiu  amjLjMuui 


‘■r? 

«*»  «*'  ^  W 

!  i*  4  4 

*  Iz  p(k)  dz  *  ave  z 

*  *  S 

or  about;  the  mean  or  average  z  or  p^  (this  notation  is  used  as 

t 

well  as  u^)  of  the  given  distribution 

n  _  2  '  o  _  2 

d2  =  J  (z  -  z)  p(z)  dz  *  -  Pj  *  ave  (z  -  zj 

M-3  =  J(z  -  z)3  p(z)  dz  =  p^  -  3Ui,M-2  +  -  ave  (z  -  z-)3 


M-4  =  J* (z  -  z)4  p(z)  dz  =  p^  -  4pxp'  +  6^i2^2  ”  3m*i4  =  ave(z-z)4 

J 

•  *  • 

While  the  dimensions  of  are  the  same  as  those  of  z,  the  dimen¬ 
sions  of  Pg,  p^,  ...  are  the  square,  cube,  ...  of  those  dimensions. 
It  is  often  convenient  to  make  this  clear  by  introducing  another 
quantity  of  the  dimensions  of  z,  together  with  suitable  dimension¬ 
less  coefficients.  This  _s  most  easily  done,  as  has  been  done  for 
a  long  time,  by  introducing  the  root -mean-square  deviation,  now 
commonly  callr  '  the  standard  deviation,  o,  where  a2  =  p2,  and  ex¬ 
pressing  \iy  p^,  ...  as  suitable  multiplies  of  a3,  a4,  .  .  .  . 

Various  systematic  notations  have  been  proposed  for  the  dimension¬ 
less  multipliers,  but  their  use  seems  to  lead  to  complexity  of  no¬ 
tation  whenever  several  variables  are  involved.  Consequently  we 
shall  here  use  the  following  rather  unsystematic  notation 
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2 

a  -  var  z 


ave 


t  “\2 

(z-z) 


p(z)  dz, 


J_  *■ 


islti 


^  *.  *  -  .  -  -<  ,  r  - 


*#  t  “Wfe  *=$1 


^  *•  a  VC  {2-?} 3  *  f(z-z}:  p  ( 2 }  cl  2 


raH  *  ave  (z-z)u  *  J(z-z)4  p(z)  dz , 


Go5  *  ave  (z-z)5  =  pz-z)5  p(z)  dz. 


We  call  7,  r,  and  G  relative  third,  fourth  and  fifth  moments. 

3.  propagation  into  Averages  and  Variance. 

If  now  we  write  yv  and  G1  for  the  correspond¬ 
ing  quantities  associated  with  w^>  and,  similarly  c^-*  J'21  ^2 
and  for  those  associated  with  Wg,  and  so  on  and  on,  we  can 
write  down  the  generalized  propagation  of  error  formulas  for  the 
case  where 

z  —  h(w^,  w 2 j  ) 

and  the  fluctuations  in  w1 ,  Wg,  w^.  are  Independent.  These 

formulas  involve  the  partial  derivatives  of  h(wx,  Wg,  wk), 

evaluated  at  the  point  where  each  w&,  a.  =  1,  2,  .  ..,  k  takes  its 
average  value.  We  denote  such  numerical  values  of  partial  de¬ 
rivatives  --  such  derivative  values,  as  we  shall  say  later  —  by 

h  h  h  hr,  h  ...  where  the  number  of  subscripts 

nl’  23J  111"  V  aab" 

shows  the  number  of  differentiations  and  the  particular  subscripts 
specify  directly  the  subscripts  of  the  w's  with  respect  to  which 
these  differentiations  were  carried  out. 

If  we  write  z  for  the  value  of  z  when  all  w  are  at 

their  average  values,  then 


■m 


~:W 

ft 


Bk 


*  * 


-  I  -*J  - 


ave  is 


A  2  h 


aa  aa 


*  B  z  haan  4  V 


+  54  2  haaaa  Va*  +  T  E*  haabb  4  V 


+ 


T50  ^aaaaa  uaV  +  I?  £*  haaabb  4V’°b 


G„  cr  ^  ^ 


3„  2 


+  terms  of  order  >  a° 

where  the  starred  summation  signs  are  to  be  interpreted  as  in¬ 
cluding  each  distinct  term  which  does  not  identify  subscripts 
once  and  only  once.  Thus  would  Include  ^2o2  once 

(not  as  ct42c752  +  o2  a2)  but  not  c44,  while  Z*o  3a  2  would  in- 
elude  a4  a 5  +  a 5  cr4  (and  not  merely  one  of  these)  but  not  a45. 
This  same  convention  on  starred  summation  signs  will  be  followed 
throughout.  (Note  that  haabb  is  symmetric  in  a  and  b  by  its 
definition. ) 

Similarly,  the  variance  is  given  by 


var  z 


^  2  2 
Z  n  a 

a  a 


+  Z  h  h  y  a  • 

a  aa .  'a  a 


(formula  continues) 


*4-2  r.  a„  ~  t-  Z  haa^(rc-l  )aah 


3 


asa  a  a 


*  z*(haftatb  +  hab2  +  haabV°aV 


4-  y-'n^  2  h  h  {Ju  ”  +  i  r  h  h  ( Q  ~  *y  \  f-  5 

12  H  SEES  EE  O  c tcL  &&&  '  J&  ? &  *  ^  \ 


■  Z*(?  hahaabb  +  I  haahabb  +  habhaab 


3  kaabhb^aaa  ab 


+  terms  of  order  >  a 
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The  first  term  on  the  right  hand  side  is,  of  course,  the  form 
which  the  classical  propagation  of  error  formulas  take  for  the 
ease  of  independent  w's. 

If  we  should  be  content  to  assume  that  the  distribu¬ 
tion  of  z  were  normal  (*  Gaussian),  then  we  might  answer  our 
target  question  easily.  For  the  tolerance  6  corresponds  to  a 
standardized  deviate  of 

(z  +  6)  -  ave  z 


v'var  z 

which  Is,  in  ease  of  normality,  easily  converted  into  the  prob¬ 
ability  of  exceeding  5  by  making  one  reference  to  a  table  of  the 
cumulative  normal  distribution.  However,  such  an  assumption 

cculd  sometimes  prove  most  dangerous,  since  for  example,  h  (w 

1  1 


w 


2J 


wk)  c.nd  hg(w^,  Wg,  ...,  wk),  where  hg  -  17(h1-zQ)^  cannot 


s 

I 


if 

1 

I 


1 

§ 


1-7  - 


both  be  nearly  normally  distributed.  We  must  be  prepared  to 
carry  our  formulas  further,  either  so  that  we  may  show  the 
normal  approximation  satisfactory,  or  so  that  we  may  do  better, 
as  circumstances  dictate.  It  Is  natural  to  attempt  this  by  using 
the  higher  moments  of  z. 

4.  Semlnvarlants  or  cumulants 

There  is  much  to  gain  by  avoiding  a  frontal  attack 
on  these  higher  moments.  First  Thiele,  and  then  R.  A.  Fisher, 
have  vigorously  pointed  out  the  advantages,  particularly  when 
effects  are  being  combined,  of  replacing  the  higher  moments  about 
the  mean  by  another  set  of  quantities,  called  semlnvarlants  or 
cumulants,  for  which  the  notation  k^,  *2,  ls  usuWl.  The 

first  two  of  these  are  the  "average"  and  the  "variance".  There 
has  been  no  general  agreement  on  similar  names  for  the  succeeding 
quantities,  but  the  word  "skewness"  has  long  been  associated  with 

Although,  unfortunately,  "skewness"  has  also  been  associated  \ 
with  various  dimensionless  quantities,  we  shall  here  call  it- 
self  the  "skewness",  with  the  natural  3-letter  abbreviation  ske  z 
It  is  the  last  cumulant  to  reduce  to  a  moment  about  the  mean 
(it^  *  M-g).  In  qualitative  terms  It  usually  tells  us  about  a  dis¬ 
symmetry  of  the  distribution  concerned  —  positive  skew  corres¬ 
ponding  to  a  "longer"  tail  --  one  more  slowly  declining  to  zero  — 
on  the  ri  t  hand  of  the  distribution  than  on  the  left  —  negative 
skew  corresponding  to  the  mirror  image  of  this  situation. 


i 

i 


* 


The  next,  fourth,  cumulant ,  *4  -  has  long 

been  connected  with  the  word  "kurtosls"  —  a  word  which  led  even 
eminent  statisticians  to  mnemonic  devices.  We  shall  use,  as  an 
interim  measure,  subject  to  a  better  suggestion,  the  term  elon- 
gatlon"  and  the  abbreviation  elo  z.  Positive  elongation  usually 
means  that  the  tails  (or  possibly  one  dominant  tall)  of  the  dis¬ 


tribution  fall  toward  zero  more  slowly  than  is  the  case  for 

(comparable)  normal  (Gaussian  or  Maxwellian)  distributions. 

Thus  distributions  with  tails  decreasing  like  e  ",  rather  tnan 
2 

like  e~cz  ,  are  almost  sure  to  have  positive  elongation.  If 
tne  tails  are  shorter  than  those  of  normal  distributions,  as  no¬ 
tably  in  rectangular  (sometimes  called  uniform/  distributions, 
or  in  U-shaped  distributions,  the  elongation  is  usually  negative. 

Various  devices  have  been  psoposed  to  make  use  of  in¬ 


formation  from  higher  moments  or  cumulants  in  selecting  better 
approximate  distributions  than  the  normal.  Three  broad  classes 
deserve  mention  here,  although  none  of  them  is  available  ade¬ 


quately  packaged  for  easy  use. 

Methods  based  on  transformation  into  terms  of  increased 
normality  were  strongly  urged  by  Edgeworth  in  the  early  years  of 
this  century.  However,  no  generally  useful  methods  were  developed. 
Transformation  In  the  special  situation  with  which  we  are  now 
dealing  is  investigated  in  Parts  IV,  VIII  and  IX,  with  mildly  dis¬ 


couraging  results . 


Methods  baaed  on  series  expansion  In  terms  of  succes¬ 
sive  derivatives  of  the  normal  distribution  have  a  respectable 
antiquity,  were  modified  by  Edge-worth,  arid  were  brought  Into 
their  most  useful  and  easily  usable  form  to  date  by  Cornish  and 
Fisher*  in  1937. 

Methods  based  on  a  family  of  curve  shapes  based  on  a 
certain  .differential,  equation  with  a  number  of  adjustable  para¬ 
meters  were  introduced  about  the  turn  of  the  century  by  Karl 
Pearson.  Practical  use  of  these  ''Pearson  curves"  has  tended  to 
be  onerous,  until  the  appearance  in  1951  and,  in  Improved  form, 
in  1954  of  convenient  tables  of  $  points**  In  standard  measure 
for  distributions  following  Pearson  curves.  In  instances  where 
a  standardized  deviate  falls  somewhere  near  a  tabulated  value, 
those  tables  offer  the  handiest  solution  so  for  available,  and 
their  use,  together  with  the  use  of  Cornish-FIsher  technique 
might  well  have  been  described  in  the  present  memorandum.  In 
view  of  the  fact,  however,  that  a  still  more  convenient  packaging 
seems  possible,  all  detailed  discussion  Is  being  omitted. 

*E.  A.  Cornish  and  R,  A.  Fisher  "Moments  and  cumulants  in  the 
specification  of  distributions"  4  Revue  de  Instltut  Intern, 
de  Statistique  1-14  (1937)  reprlnTed  as  paper  l*i  EfT  A. 

sher.  Contributions  to  Mathematical  Statistics,  New  York, 
Wiley,  1$5^ 

**Maxine  Merrington  and  E.  S.  Pearson  "Tables  of  the  5$  $nd 
0.5$  points  of  Pearson  curves  (with  argument  f>i  ard  62)  ex¬ 
plored  In  Standard  Measure"  38  Blometrika  4-10  (1951)  also 
In  Improved  form  as  Table  42  of  'Biometrl¥a  Tables  for  Statls- 
t 1 clans ,  Volume  1,  E.  S.  Pearson  and  H.  0.  Hartley  editors, 
Cambridge  University  Press  195^. 


Either  the  Cornish -Fisher  or  Herring ton -Pears on  tech¬ 
nique  starts  conveniently  from  dimensionless  ratios  of  cumulants. 
Those  entering  into  the  Cornish-Pisher  formulas  are 


and 


?1 


ske  z 


'var  z  ] 


37? 


3 

17? 


^2 


elc  z  _  ^4 


both  of  which  vanish  for  the  normal  distribution.  Those  enter¬ 
ing  into  the  Me rrlng ton- Pearson  tables  are 


[ske  z  ]2  =  2 

[var  z  ? 1 


44 

- *5  ~  0 

*2 


-f 


elo  z 

"  - - n- 

[var  z 


3  +  72* 


whose  values  for  a  normal  distribution  are  0  and  3.  Both  ap¬ 
proaches  make  use  of  the  standardized  deviate,  or  deviate  ex¬ 
pressed  In  standard  measure 

(zQ  +  6)  -  (ave  z) 

/var  z 

relating  this  either  to  a  modified  deviate  which  may  be  referred 
to  a  unit  normal  distribution,  or  directly  to  the  probability  of 
exceeding  z  +  6. 


-  X-il 


Te  keep  the  record  cAear,  we  must  emphasize  that  both 
methods  are  only  approximate,  though  in  practice  usually  ade¬ 
quate.  Experience  with  the  Pearson  curves  near  the  $  points  so 
far  tabulated  has  been  excellent.  In  further  view  of  the  con¬ 
venience  of  the  Merrlngton-Pearson  tables,  their  use  for  true 
probabilities  of  exceeding  tolerance  between,  say,, 0.3$  and 
$$  is  at  present  to  be  recommended.  For  probabilities 
greater  than  about  10$  (and  leas  than  90$)  Comish-Fisher  treat¬ 
ment  may  be  r*at.  In  either  case,  at  the  user's  choice,,  the 
purpose  of  more  detailed  calculation  may  be  either  to  indicate 
the  adequacy  of  the  normal  approximation,  or  to  provide  a  better, 
satisfactory,  approximation. 

5.  Propagation  into  Skewness  and  Elongation 

We  can  now  state  the  analogous  fonnulas  for  the  leading 
terms  in  ske  z  ,  and  elo  z  ,  namely 


ske  z  *  2  h  ^  7  a  3 

a  'a  a 


+  t  2  ha2haa'ra-h°a4  +  6z* 


+  V  2  ha2|WV>a)oa5  +  i  s  Vaa^V2^ 


+  3  2*  V-ab^a-^aV 


+  |  S*(ha2habb  +  2haahafchb  +  2Vaabhbha°a3ob2 


+  terms  of  order  >  a 


-  x-ia  - 


and 

elo  2  ~  Z  ha^ra“^ffa^ 


2  *  ha3haa<V4  ’’aK5 


+  terms  of  order  >  a 


The  detailed  correctness  of  the  results  of  this  ac¬ 
count,  both  above  and  in  sections  to  ccme,  owes  much  to  careful 
checking  and  error  finding  by  Miss  M.  S.  Harold.  Thanks  go  to 
her  from  the  writer,  both  on  his  own  behalf  and  o i  behalf  Ox 

future  users. 
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TSaMS  FOR  INDIVIDUAL  V’ARTARr.F'.c! 


6.  General 

We  come  now  to  the  uses  c.f  changes  in  the  terms  of  study. 

The  fraction  of  some  population  of  electronic  assemblies  with  a 

certain  frequency  greater  than  400  cycles  per  second  is,  of.courae, 

exactly  the  same  as  the  fractions  (1)  with  the  square  root  of  the 

frequency  greater  than  20  ( cycles/sec J1/2,  or  (11)  with  the  common 

logarithm  of  the  frequency  (In  cycles  per  second)  greater  than 

2.60206-.  If  it  IS  easier  to  work  with  one  of  these  equal  fractions 

rather  than  either  of  the  others.  It  Is  clearly  to  our  advantage  to 
do  so. 

In  some  circumstances,  where  the  use  of  unfamiliar  terms 
of  analysis  Is  called  "transformation",  or  even  "transforming  the 
data'  there  13,  to  the  minds  of  some,  an  unfortunate  flavor  of 
cooking".  While  that  feeling  Is  almost  always  (If  not  always!) 
unjustified,  in  the  present  situation  there  is  no  slightest  excuse 
for  a  similar  feeling.  Suppose  that  Individual  assemblies  are 
nested  --  for  attenuation,  resistance,  amplification,  critical  fre¬ 
quency,  or  what  have  you.  And  suppose  further  that  the  result  is 
displayed  by  the  motion  of  a  needle  across  a  scale.  Certain  points 
on  that  scale  correspond  to  assemblies  below  the  tolerance,  others 
to  assemblies  above  tolerance.  And  this  need  be  changed  in  no  way 
whatsoever  If  the  numbers  on  the  scale  represent,  for  example: 

(1)  ohms 
(ii)  y/oEmB 


(ill)  logarithm  of  ohms 
(iv)  (ohms}"'1/2 
(v)  (ohms)"1  »  mhos 

or  any  other  similar  scale*  Juet  so  long  as  the  needle's  position 
exactly  at  the  tolerance  is  Kept  the  same.  Just  so  long  will  we  be 
dealing  with  the  same  problem. 

All  tnis  has  referred  to  changes  in  terms  from  z  to  some 
function  y  of  z,  The  situation  in  a  change  of  terms  from  some  w 

a 

to  some  function  v&  of  wa  Is  similar,  but  not  identical.  If  we 
are  to  use  either  the  classical,  or  the  generalized,  propagation 
formulas,  we  need  to  know  something  about  the  moments  of  the  dis¬ 
tribution  of  w&  in  one  instance,  and  those  of  the  distribution  of 
va  on  the  other.  This  fact  seems  to  have  worried  persons  with  a 
mathematical  background  and  orientation  far  more  than  it  should  -- 
as  is  rather  natural.  If  one  starts  to  treat  a  situation  by  say¬ 
ing  Let  us  assume  that  we  know  the  low  moments  of  the  w  "  and 

a  1 

presume  that  "we  may  treat  the  Joint  distribution  of  the  w  as 

a 

multivariate  normal",  and  sticks  rigidly  to  these  hypotheses,  as 
would  be  appropriate  in  pure  mathematics,  then  he  must  realize  that 
he  does  not  know  the  low  moments  of  the  v  exactly,  and,  moreover, 
if  the  Joint  distribution  of  the  w^  was  exactly  multivariate  normal, 
then  that  of  the  v^  cannot  be .  If  the  first  set  of  starting  hypoth¬ 
eses  was  taken  as  gospel  engraved  on  tablets  of  stone,  then  the 

| 

transformation  appears  to  weaken  and  distort  the  availablfe  knowledge. 
But  in  a  practical  situation  this  is  only  an  appearance,  One  never 
knows  the  low  moments  exactly,  and,  although  one  may  ha$T  fairly 


goo  J  liEsteSjr  one  often  doe®  not,  (In  tolerance  analysis  situa¬ 
tions  the  low  moments  are  often  of  a  "but  what  if"  character,  and 
are  far  from  precise.)  Moreover,  if  approximate  normality  could  be 
assumed  for  the  w ‘ s,  It  is  not  unlikely  that  approximate  normality 
can  be  equally  well  assumed  for  the  v's.  We  must  never  forget 
that,  in  applications  of  mathematics,  exact  hypotheses  are  usually 
only  approximations  to  the  real  situations  and  that  other  approxi¬ 
mations  may  be  as  good,  or  even  better,  as  those  we  first  made. 

This  apparently  subversive,  but  actually  usually  unim¬ 
portant,  effect  of  transformation  of  Individual  variables  Is  enhanced 
in  dealing  with  propagation  fomrulas .  These  formulas  Involve  de¬ 
rivative  values  at  "the  average  point",  and  the  average  value  of 
va,  where  va  is  a  given  function  of  w&,  is  not  exactly  the  same  as 
the  given  function  of  the  average  value  of  wg .  Thus  the  average 
point  shifts  under  the  transformation  and  the  derivative  values  we 
need  are  not  only  for  derivatives  with  respect  to  new  variables, 
which  we  can  evaluate  by  simple  transformations,  but  are  evalua¬ 
ted  at  a  new  point,  if  we  had  known  the  old  average  point  exactly, 
then  this  loss  In  firmness  would  have  been  relevant,  though  prob¬ 
ably  not  serious ,  But  we  usually  do  not  know  it  exactly,  and  In 
practice  we  are  not  appreciably  worse  off. 

Thus,  while  it  Is  conceivable  that  we  might  be  able  to 

handle  the  moments  of  the  w  and  evaluate  derivatives  at  the  w  - 

a 

average  point  --  and  still  not  be  able  to  do  the  same  in  terms  of 
the  va  —  such  a  situation  is  most  unlikely  to  occur.  In  practical 
situations,  transformation  of  individual  variables  may  usually  be 
undertaken  quite  freely. 


II-4 


Trans  format 


-9/  -th&  Individual  ( Component  j  Variables 
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We  can  discover  these  transf ormatlons  by  complex  argument* 

If  we  wish,  merely  by  letting  w,  be  some  function  of  v  ,  and  find- 

^  a 

ing  successive  conditions  on  the  derivatives  of  this  function  eo 
make  more  and  more  terms  vanish.  nd  then  finding  a  function  with 
these  derivatives.  But  we  may  also  avoid  all  this. 

Consider  what  happens  if  all  w's  except  w  are  at  their 

a 

average  values,  then  s  depends  on  w  alone,  but  usually  not  in  a 

a 

linear  way.  If  we  may  reasonably  define  v„  so  as  to  make  this  Hp- 

.  . "■  -  a. 

pendence  linear,  then  in 

z  =  s(v1,v2J  . . .,vk) 

the  higher  urenixed  derivative  values  gaa,  eaaa,  Saaaa ,  which  we 
shall  later  refer  to  as  successive  derivative  values,  will  all 

vanish,  and  the  formulas  will  be  simplified. 

n 

When  may  we  reasonably  dc  this?  Anytime  that  z  is  a 
|  strictly  monotonic  (ever-increasing,  or  ever-decreasing)  function 

Iof  each  w&  over  its  normal  range  of  variation.  For  if  we  write 
VV  w2' ••■*wk  for  the  average  values  of  Wg, ...,wk  and 

*a("a)  -  wa,  wa+b,...,  wk) 

I 

for  the  value  of  z  =  g(w1,w2, . . . ,wk)  with  all  the  variables,  other 
than  wa,  at  their  average  values,  then  ^a(w&)  is  naturally  described 

ias,  and  called,  the  partial  effect  of  w&  at  the  average  point  (really 

along  the  line  through  the  average  point  for  which  w., ,  w„3 ,  w 

l'  2  * "a-1* 


|w>!,*WMir"||i|iw*,iiwia,>i,tc^^ 


*  msmrnm: 


■  Z\  . 

.  -i.  . 


wa+l'/  *  *  ‘ ' wk  w  1  s  except  w&  --  are  constant ) .  if  we  could 

Introduce  va(lj  *  ^{w^}  as  a  new  variable,  so  that 


z 


a) 


describee  the  response,  then  the  successive  derivatives  (after  the 
first)  of  g(1)(-,-,...,-)  with  respect  to  vj1'',  evaluated  at  the 
v-values  corresponding  to  the  w-average  point  will  vanish.  This 
will  be  so  because  z  =  vj1)  so  long  as  the  v^'s  other  than  vj1) 
are  at  the  values  corresponding  to  this  point,  since  the  w's  other 
than  wa  will  all  be  at  their  average  value.  Two  difficulties  then 
remain:  First,  the  transformation  may  not  be  permissible.  Second, 
we  wish  the  successive  derivatives  to  vanish  at  the  v-average  point 
rather  than  at  the  w-average  point. 

The  transformation  will  fail  to  be  permissible  if  two 
values  of  „a  correspond  to  the  same  value  of  ^(wj.  In  practlce 

we  require  such  single-valuedness  of  the  Inverse  only  for  varies 
of  wa  reasonably  near  wa>  Such  a  failure  of  uniqueness  will 
ordinarily  only  occur  when,  because  of  design  (or  perhaps  by  acci¬ 
dent)  the  nominal  value  of  w&  has  been  chosen  to  make  z  exactly 
or  nearly  a  maximum  or  minimum  with  respect  to  changes  in  w 
In  most  such  situations  ha  will  be  small  and  h^  will  be  important. 

It  will  then  be  possible  to  choose  v&  so  that,  while  g  does  not 

aa 

vanish,  gaaa  =  6aaaa  =  ...  »  0.  (Details  are  left  to  the  reader.) 

With  this  exception,  then,  the  transformation  will  ordi¬ 
narily  be  permissible,  and  we  need  only  be  concerned  that  the  v(l>- 
average  point  Is  not  the  same  as  the  *r -average  point.  If  we  did  not 


li~6  - 


know  Che  location  of  the  w-average  point  precisely,  this  fact  will 
not  concern  us  very  much.  Either  we  did  not  know  the  partial  effect 

very  well  either,  since  it  depends  on  the  average  point,  and  ever^- 

thing  la  really  about  eauallv  ,,  .  . 

J  uuo  equally  ha*.y,  or  we  were  able  to  work  with 

in  whici  the  coordinates  of  the  w-average  point  is  ex¬ 
pressed  by  letters  which  we  may  substitute  for  when  we  are  done  -- 
and  nothing  prevents  us  from  making  a  substitution  corresponding 
to  the  new  average  point,  although  (if  „e  have  exact  knowledge) 
we  may  have  to  Iterate  a  little  to  find  the  correct  values  for  sub¬ 
stitution.  Thus  the  shift -in-average -point -problem  is  not  likely 
to  be  important. 

If  we  follow  the  program  suggested  above  in  detail,  we 
would  reach 


with  not  only 


z  “  v2,  . .  .,v,  ) 


®aa  ^aaa  • • •  ~  ® 


but  also 


ga  »  1 


for  all  a.  it  is  often  natural  to  introduce  particular  Vs  which 
make  the  successive  derivatives  vanish,  but  do  not  have  all  g  =  1. 
W e  teSln*  then,  with  this  more  general  case. 

— Reduced  Propagation  Formulae 

In  the  usual  Instances,  we  are  merely,  in  system-component 

terns,  planning  to  measure  the  performance  of  a  component  by  the 

P  ormance  of  a  system  containing  this  component  but  with  all  the 

other  components  selected  to  be  at  the  mean  values  of  each's  dls- 
trlouuion. 


*iik  ikii  i?slli 


¥ 


'I 
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Uider  these  conditions,  the  formulas  simplify  consider¬ 
ably.  We  new  have  z  =  vgJ  .  ..,vk),  where 

(1)  the  V(  fluctuate  independently, 

(2)  gaa,  gaaa ,  etc.  all  vanish, 

and  the  generalized  propagation  formulas  reduce  to  those  of  Table  1. 

Notice  that,  although  the  first  correction  terms  to  the  average, 

variance  and  skewness  are  all  of  order  o^,  their  relative  orders 
4  2 

are  a  ,  a  and  a,  respectively.  Thus  correction  terns  are  most 

likely  to  have  nade  a  noticeable  fractional  change  in  the  skewness, 

next  in  the  variance,  and  least  in  the  average.  Note  further,  that, 

except  in  the  leading  term  for  the  elongation,  only  the  variances 

aa  and  skewnesses  70^  of  the  individual  variables,  the  v  ,  enter 

a 

the  generalized  propagation  formulas  in  terms  of  order  <  0^.  Thus 
the  requirements  as  to  knowledge  of  component  variability  (in  the 
He,w  terms' )  are  not  as  stringent  as  we  might  have  feared. 


Table  1 


Propagation  in  the  Special  Case  Where  g„.  etc .  Vanish 

3.3.  3.33 

(Individual  Variables  measured  proportionately  to  partial  effects) 


ave  z  =  z. 


(i.e.  z  calculated  at  the  average  point) 


+  V  **  Saabb  °a  V 

1  v*  „  3  2 

H?  ®aaabb  ya0a  ab 


+  terms  of  order  >  a 


var  z  = 


„  2  2 
2  sa  aa 


(i.e.  the  classical  term) 

2  _  2 


+  2*  ^sasabb  +  ®ab  +  saab®b^aa  ab 
+  2*  ^  ®a®aabb  +  sabsaab  +  "J  ^aaabs,b^aaa  ab 


+  terms  of  order  >  a 


ske  z  = 


2  sa3 

+  6s*  Wab°aV 
+  32*  Mab^a^WV 

+  |  2*  (sa2Safcb  +  2SaS3.abgb)’'a0a3ab2 


4-  terms  of  order  >  a 


elo  z  = 


2g  ^(T  -3)  a  21 
v  a  a 


+  terms  of  order  >  a 


o 


(Symmetry  of  g  .  and  g  ,,  in  a  and  b  to  be  recognized  in  interpret¬ 
ing  2*) 
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We  are  now  prepared  to  -seek  further  i  ondensatlon  of  these 
propagation  formulas.  We  may  gain  a  little  by  arranging  for  all 
the  ga  to  be  unity  --  we  had  already  expressed  the  individual  vari¬ 
ables  in  terms  such  that  gOQ,  etc,,  all  vanish,  we  can  still 

choose  the  size  of  the  units  in  terms  of  which  they  are  expressed. 
We  shall  then  be  measuring  each  individual  variable  exactly  in 
terms  of  its  partial  effect  of  its  response.  (This  is  entirely 
analogous  to  Inspecting  components  by  measuring  the  performance  of 
a  standard  system  in  which  all  other  components  are  at  the  average 
values  of  their  distributions . ) 

If  we  introduce  reduced  or  standardized  quantities  t 

d. 

through 

Ta  *  Sa  aa  , 

2 

then  t  will  be  the  variance  of  the  reduced  variable  x  =  g  v  . 
a  a  ea  a 

Moreover,  y  and  F  will  be  the  same  for  the  distribution  of  x 
a  a  a 

as  for  that  of  v&.  We  only  need  to  define  some  t's  by 

1-  -  u  _  ^abbb 

ab  -  gagb  ’  abbb  -  ^7? 

and  the  like,  where  the  subscripts,  as  on  t's  and  s's  generally, 
mean  simple  differentiation,  in  order  to  obtain  the  pro¬ 
pagation  formulas  in  the  form  given  in  Table  2, 


i 

avs 


var 


ske  i 


elo  z 


(Symmetry 
ing  s#. 


me  Propagation  Formulas  in  Reduced  Perm 
Individual  variables  measured  as  partial  effects) 


(i.e;  ?.  evaluated  at  the  average  point) 


+  f?  *a.bbi*a\ 


2_  2 
v 

hJ 


2V  )tTT  t.  )  “V  <r  3.,  C. 


I?  aaabb  j  ^aTa  Tb 


+  teims  of  order  > 


z  —  St 


a  (l.e.  the  classical  term) 

2 

aab  "  uabb)  Ta  Tb 


*  \  fcab2  +  taab  +  ^  ''  2~  2 


] 


k  laabb  +  tabt'aab  +  ?  fcaaab  K^b 


3.  2 


+  terms  of  order  > 


Z  y  t  ■ 
a  a 


(i.e.  the  sum  of  component  skewnesses) 


+  2*  (6  t  ,  )T  2t  2 
ab'*a  b 


+  S*  <3  ‘ab^V^aV+^wft abb)7^a\2 


+  terms  of  order  >  t 


o 


2  (r  -3)t  ‘ 4 


(i.e,  the  sum  of  component  elongations) 
+  terms  of  order  >  T° 


of  fcab  and  taabb  ln  a  and  b  t>e  recognized  In  Interpret- 


.■i 


-  11-2 i  - 

In  preparing  Table  2  we  have  tried  to  throw  each  for¬ 
mula  into  Its  most  usable  fora,  and  have  consequently  combined 
numerical  coefficients  with  the  t!s  they  modify,  in  general,  the 
most  Important  terms,  beyond  the  leading  term  in  each  formula, 

seem  likely  to  be  the  second  and  third  terms  for  the  skewness  and 
the  second  term  for  the  variance. 

— Ifeg-.Causes  and  Justifications  of  Slmplif lcatlon 

We  have  now  reached  a  fairly  compact  set  of  propagation 
formulas  which  contain  all  teras  through  those  of  order  .  it 
is  time  for  us  to  inquire  into  the  assumptions,  explicit  and  im¬ 
plicit,  which  have  entered  into  the  development-indicating  the 
justification,  complete  or  partial,  of  each.  The  list  to  be  ex¬ 
amined  Includes  assumptions  that: 

(i)  the  functions  considered  are  sufficlentlv 
difi erentiable  to  have  convergent  Taylor 
series, 

(ii)  the  functions  considered  are,  locally, 

sufficiently  like  polynomials  to  be  reason¬ 
ably  well  approximated  by  the  first  terms 
of  their  Taylor  series, 

(iii)  the  variables  taken  as  individual  (=  component) 
variables  have  independent  fluctuations  (at 
least  the  cross -moments  of  order  <  5  vanish 
when  the  variables  are  expressed  in  particular 
terms ) , 


r  r-i  o 


I 


ilv)  we  may  use  axi  expansion  around  the  average 
values  of  the  Individual  («  component)  vari¬ 
ables,  expressed  in  whatever  terms  we  may  use, 

(v)  the  response  (-  overall)  variable  is  a  strictly 
monotone  function  of  each  Individual  (*  com¬ 
ponent)  variable. 

The  strength  of  (i)  and  (li)  ran  be  fairly  well  discussed 

in  tonns  of  z  =  h(w)  where  the  function  h  is  useu  as  a  rectifier. 

If  we  have  a  quadratic  rectifier  7  -  , 

xi. er,  z  -  cw  ,  then  the  basic  formulas 

(though  net  those  reduced  by  transformation  of  Individuals)  will 
apply.  If  „e -have  an  Ideal  linear  rectifier,  z  =  bjw|,  then  (1) 
falls  because  of  the  misbehavior  of  the  derivative  at  zero.  If 
we  observe  that  actual  "linear"  rectifiers  deviate  from  'Ideal' 
behavior  at  zero,  and  manage  to  behave  as  If  they  have  a  continuous 
derivative,  we  find  that  (i)  Is  satisfied,  but  that  we  are  likely 
to  be  troubled  with  (11),  since  the  first  derivative  Is  nearly 
constant,  except  near  zero,  where  It  changes  rapidly.  This  may 
well  be  far  from  polynomial  behavior'  Only  If  our  'linear'  rec¬ 
tifier  is  quite  far  from  'ideal'  can  we  expect  even  generalized 
propagation-of -error  formulas  to  be  a  close  approximation.  The 
application  ol  (i)  and  (11)  to  other,  more  complex,  situations 
tendsto  follow  the  pattern  set  it,  this  example.  If  a  single 
vaxue  of  each  first  partial  derivative  Is  sufficient  to  guide 
even  the  roughest  design  (and  sometimes  when  this  Is  not  true) 
oth  (1)  and  (11)  are  likely  to  be  well  enough  satisfied. 


> 
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The  classical  formulas  are  frequently  written,  not  only 
for  Independent  variables 

,ar  »  '  2  ffsH  var  xa 

a 

but  also  for  the  general  case 


+  22*[s^l  oov  (xa’V 

Why  have  we  not  followed  the  pattern  in  setting  forth  the  general¬ 
ization?  Both  because  it  would  be  difficult,  and  because  it  would 
really  add  little.  The  added  complexity,  for  the  classical  formula, 
is  small,  and  pays  for  itself  in  flexibility  —  in  the  ability  to 
work  with  Intermediate  variables  whose  variances  and  convariances 
have  been  already  obtained  by  previous  uses  of  propagation  formulas 
--  in  the  ability  to  deal  with  situations  where  it  is  hard,  or  in¬ 
convenient,  to  find  statistically  independent  variables.  In  prac¬ 
tice  it  seems  to  be  the  first  ability  rather  than  the  second  which 
Is  usually  important.  Thus  the  availability  of  the  classical  for¬ 
mulas  for  the  dependent  case  buys  something,  but  not  very  much. 
Building  provision  for  statistically  dependent  component  variables 
into  the  generalized  formulas  would  greatly  complicate  them  --  so 

much  so  that  it  seems  unlikely  that  terms  beyond  those  of  order 
3  4 

(or  perhaps  t  )  would  be  written  out.  Thus  the  effort  of  seek¬ 
ing  basic  individual  variables  with  independent  fluctuations  is 
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likely  to  be  Tar*  less  than  the  effort  of  working  with  formulas 
for  dependent  ones .  Thus  (ill)  has  much  practical  Justification. 
(In  those  instances  where  it  does  not  hold,  it  will,  of  course, 
be  necessary  either  to  use  the  classical  formulas  or  develop 
further  approximations.) 

What  of  the  remaining  assumption  --  that  we  use  an 
expansion  around  the  average  values  of  our  individual  variables? 
When  we  recognize  the  inaccuracies  of  our  knowledge  in  any  prac¬ 
tical  application,  this  assumption  does  not  really  seem  serious . 

If  our  knowledge  is  basically  mathematical  and  approximate,  gained 
by  calculating  values  of  partial  derivatives  at  a  certain  point, 
then  it  may  be  true  that  this  point  would  not  be  the  average  point 
in  a  real  situation.  But  it  is  quite  likely  to  be  true  that  we 
do  not  know  accurately  where  the  real  average  point  will  be.  This 
means  that  we  must  Interpret  cur  derivatives,  and  all  that  hangs 
upon  them,  with  caution  and  some  allowance  for  variation.  This 
will  be  the  case  for  any  formula,  whether  or  not  It  expands  around 
the  average  point. 

In  order  to  deal  wisely  with  this  possibility,  we  need 
to  know  how  the  g&  and  tab  change  when  we  move  from  the  point 

^10*  v20'  *•••’  vk0  t0  the  nearby  Point  vn,  vgl ,  .  .  . ,  vkl .  If  we 
let  v^,  Vg, ...,  vk  be  the  displacements,  as  measured  by  the  re¬ 
sponse  at  the  average  point,  i.e. 


g(v10J v20>  *  *  '*  vaP,,,,v! rf))-S( 


al' 


'kO- 


r10'v20'  *  *  *  *  va0#  *  *  *,vkC>) 


-V,, 


.  *'  -  ••  '-f 


•  -C 1 


*  XI*1$  " 


then  the  original  forms  and  first  order  correction  terms  am  In¬ 
dicated  by  the  following  expressions 

+  £c  V  vo  + 

^ab^  "r  ^'abc”tac'*^bc ^  vc  +  *  ‘  *  ^ 

"'aab^1  +  ^ toabc~^  tac  ~  tbc^vc  +  *  *  *  ^ 

and  so  on,  including 

taaabb^1  "  2c ^fcaaabbc  ”  ^  ^ac  "  ^  ^pc^vc  +  *  *  * )  * 

These  formulas  provide  the  necessary  guide,  should  it  ever  be 
necessary  to  quantitatively  judge  the  adequacy  of  the  approxima¬ 
tion  with  which  (iv)  holds.  We  may  expect  that  it  will  usually 
be  good. 


IS*#*** 


Ill  EXAMPLES 


10.  The  First  Example 

‘fie  now  report  fene  results  of  quantitative  discussion 
of  two  moderately  simple  examples  beginning  with  the  time  delay 
of  a  lumped-constant  delay  line,  each  section  of  which  consists 
of  an  inductance  L  and  capacitance  C  (initially)  of  the  same 
nominal  values,  (The  detailed  Calculations  may  he  found  in  Part 
VI.)  The  orig  nal  functional  relation  may  be  written 

delay  -  z  =  +  +  ...  W  % 


$ 

1 


1 

| 


where  we  have  used  odd  subscripts  on  inductances  and  even  ones 
on  capacitances  so  that  a  subscript-  points  uniquely  to  a  com¬ 
ponent  , 

For  convenient  Insight,  we  fix  the  units  for  L  and  C 

so  that  their  nominal  values  are  both  10,  This  means  that  a  =  1 

a 

corresponds  to  a  standard  deviation  of  10#  In  the  inductance  or 
capacity  of  component  a,  (Realistically,  then,  o!s  may  be  as 
large  as  2,  but  are  more  likely  to  be  fractions  like  0.6,  0.2  or 
0 . 06 . ) 

Direct  application  of  the  original  formulas  now  gives 
(terms  arranged  In  order  of  size  of  coefficient  --  those  with 
coefficients  smaller  than  0,000005  omitted) 


ave  z 


lOj  -  0.0125  £  a&  +  0,00062  2  7^0^ 


4 


-0.00002  ^(2r21_1cr2i_1  +  2r2ia2i4^ 


+  ,  * 


+  0.00016  X 


(5r!i-r1)a21.14  ..  40 4  (5r81-i)o21‘ 


4  0.00002  2(0a-^)aa  -O.OOOuJ  Z(?sl_102i.1+721as,)c2 


2r,  2 

21/u2i-l  a21 


+ 

'  »  * 


ske  z  =  0.125  Z  y  a 

'3.  3 


-  0.00937  I. 


^r2i-l_1  ^c2i-l  ‘  4a2i-l2  °2i2+^r2i_;i^0: 


21  '  21  '  u2i 


+  0.00023  2  (3G„-47  )a 

c.  a.  d 


-  0.00047  2 


(3^i-l+2)°2i-l  +  (3721+2)cr21 


2  2 

a21-l  a21 


4* 

•  *  *  $ 


elo  z  =  0.0625  S  (ra-3)ca4  +  0.00625  2  ((^-4^)0  5 


+  ... 


and  It  Is  clear  that.  In  every  case,  at  most  the  first  correc¬ 
tion  term  will  have  any  numerical  Importance. 

If  we  transform  the  individual  variables,  to  make 
partial  effects  linear  of  coefficient  one,  then  we  have 


where  aa  is  the  standard  deviation  of 


JU.1L. 

the  new  a  1  variable,  and 


where  subscripts  a  run  from  1  to  2j  and  subscripts  i  from  1  to  j. 
The  qualitative  simplification  is  obvious.  One  term  instead  of 
5  plus  higher  terms  in  the  average.  Two  Instead  of  five  plus 
higher  In  the  variance.  Three  plus  two  higher  instead  of  four 
plus  higher  In  the  skewness.  One  plus  higher  instead  of  two  plus 
higher  in  the  elongation.  The  quantitative  simplification  in  the 
variance  is  also  important.  If  o  <  1  for  all  a,  the  second  term 
In  the  variance  is  less  than  one -thousandth  of  the  first. 

In  the  skewness,  however,  the  second  term  retains  con¬ 
siderable  Importance  In  comparison  with  the  first.  Practical 
use  of  these  formulas  would  probably  Involve  one  tern  each  for 
average,  variance,  and,  if  it  were  required,  elongation.  If  the 
skewness  were  required,  at  least  the  second  term  would,  and  the 
third  term  might,  need  examination  In  terns  of  transformed  indi¬ 
vidual  variables,  then,  terms  beyond  classical  propagation  of  error 
need  be  considered  in  at  most  one  formula. 


average  value  cf  V^2\  -1 


standard  deviation  of  \Z'C2^  t ^ 

"p  _  .  —  ^ 

average  value  of  i 


In  terms  of  these  quantities,  the  formulas  become 


ave  z  =  Zt±3 


var  z  =  St12(Tl21.12  +  T]2JL2  +  T^2  T]^2), 


ske  z  *  72i„1Tl21.13  +  >2iM2i3  +  UTl2i-l  ^21 


+  3^2i-l"1^Tl2i-l3T1212  +  ^72i~1)T12i3T121-l  ^ 


4-  terms  of  order  >  t  , 


elo  z  =  Ztt  (r2i_i-3)Tl21_n4  +  ^r2i"3)Tl2iif 


+  terms  of  order  >  t  , 
and  are  clearly  easily  manageable , 


la,  me  second  example. 

The  attenuation  caused  by  a  ir-sectlon  network  of 
series  resistance  Rg  and  shunt  resistance  and  R,,,  operating 
between  image  impedances  R,  is 


z 


1  + 


— -  + 

2F  £ 


1 

2 


*3 


S  i 

i  ( 

R 

4*  4r 

Li..,.  1  +  \  1 

>R"  " 
\  ** 

M 

-r 

1ra  rc  rarc) 

If  we  write  this  in  terms  of  changes  from  ncminal,  the  unit  for 
each  change  being  10 %  of  the  corresponding  nominal  value,  we  find 

ra  =  r3-tt  d  +  “-1  V 


( 1  4-  0.1  Wg ) 


RC -  R  m  d  +  °-1  «c> 


where  a  is  the  design  attenuation,  and  calculate  the  necessary 
derivatives,  the  generalised  propagation  formulas  become,  if  we 
introduce  abbreviations  0,  y  \ not  the  same  as  y c  for  the 

w's)  and  5  by 


a  _  a  a  -  1 

p  2  a  +  1  ’ 


( a  -  1 ) 

2 (a  +  1 } 


the  following: 


«.  I  ~  6 


►  *  m. 
*  *  V 


*  -  a  iproviae*?  th*»  of  w 


*  w„  vanish) 


+  0.01  e(aA'  +  ff  2) 


0.001  a(yA  <JA3  +  >coc-) 


+  0.0O01  P(r(‘  +  r>„4) 


A  A  C  C 


+  0.0001  yOp^OQ^ 


-  0.00001  3(GAaA5  +  G^5) 


+  0,00001  y  (yAaA  +  7cac)  o^ac? 


+  .  , 


2  f  /-r  ^  l  ;i  _  2  .  2  * 


var  z  =  0.01  3“  (a/  +  4a/  +  a/) 


-  0,002  £2  (rAaA3  +  7cac3) 


+  0,0001  3‘ 


(3Fa  -  ^  aA4  +  (3rn  -  1) 


A  WA 


+  0.0001 


(6^  +  436 )(aA2  +  a  2)aR2+(72+437)G'  20p2 


A  '  /UB 


A  C 


-  0.00002  3‘ 


l2GA  '  7A^aA5  +  (2Gn  "  7p) 


C  'c'  UC 


(formula  continues) 


-lyryT  - 


0.00001  +  4p*)(yAoA3  +  rcac3)GB2 


-(2y  +>)Py)fj'AaA+7coc}  aACac? 


3ke  2 


0.001  ^  <^3  .  8^3  +  3) 


-  0.0003  _(rA  -  1)  0/+  (r  .  1}  4 


+  0.0012  ec6  (akd  +  Cc?)  o^2  +  0.0006  fs2y  0,Sa.  2 


A  UC 


0.00003  e3  |_(2Ga  -  3 yk)  a*  +  (2G0  -  3,c)  o/j 


0„ 00003  36^  (y  -1)^3?  q /  v  3/  2  2x 

_  7A  °A  aB  ‘  2^B  *  J)aB  (°A  ‘°C  ) 


+  (7 


C  ‘  ^cVl 


0.00003  Sr2  [{?A  -  1)„A  +  .  1)oJ  0  20  2 


0,00006  3^6  (4>  a  ^  /r  2  0  S  2  2 

!_  Va  +  4j'c°c  )ob  2^b°b  ^°a +0c 


-  0.00015  (yAaA  +  ycV  ,£  2 


iMsa .  A*  23M&SL 


v,-'  •:  ■  #ra.  \k4, 


•'  V  -  ■*  * 


^3-r ' .t  ‘4jTp 


-1 
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eio  z 


fS, 


O.OOOl  ft4 


i)ah  *  i6(rB  -3)  s.u  *  (rc-3)cc;< 


-  0.00004  g4 


4’'a)  "a5  +  <3C 


^7C) 


4*  ,  *  . 


lvhlle  some  of  these  formulas  appear  rather  lengthy,  it  is  clear 
that  the  first  two  terms  will  give  either  the  average  or  the 
variance  to  working  accuracy. 

If  we  seek  to  transform  the  individual  variables,  we 
are  Immediately  led  to  begin  by  replacing  the  resistances  R  of 

the  shunt  arms  by  their  conductances  Y.  The  attenuation  then 
becomes 

z  '  1  +  m  Kb  +  I  ^  K  +  YCj  +  f  K  +  Yc  +  Wcj 

and  we  place 


V  =  i  Q  “•  1 

A  R  o^TT 


R 


a2  -  1 


"5H- 

v  _  1  Q  ~  1 

C  R  aTT 


■  (1  +  0,1  vA) 

(1  +  0,1  vB) 

(1  +  0.1  vc) 


(Note  that  dW^/dV^ )-(dW^/dV!-, )  -  -1,  not  +1.  This  deviation  from 
our  general  practice  Is  a  matter  of  convenience  only.) 

The  generalized  propagation  formulas  now  become: 


V 


X-Jv&stom'X.-C 


*  I®  •  a  isife 


*  iix-9  - 

ave  2  *  a  +  ...  (provided  the  average  values  of  vA*  vg 

and  v  vsr  ten) 

c 

var  z  =0,01  $2  ( o ^  +  4og2  +  tfc2) 

+  0,0001  5 2(ga2  +  ac2)  aB2 
+  0.0001  72  aA2  ac2 

+  «  «  • 

ske  z  =  0.001  03  ( 7^\aA^  +  3?BaB^  +  ^CcC^ 

i-  0.0012  p25  (aA2  +  ac2)  afi2  +  0.0006  p27  oA2ac2 

+  0.00003  062  (yA  -  l)cA3  +  ( 7C  -  1)  ac3  Cg2 

+  0.C0003  P?2  (7a  -  1)<JA  +  (7C  -  i)^C  aA2aC2 

+  0.00006  062  ( 7B  -  1)  a33  ( a A2  +  ac2) 

+  ... 

elo  z  =  0.0001  (34  (rA  -  3)oa4  +  l6(rB  -  3)ob4  4-  (rc-3)ac4 

+  ... 
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mile  the  second  terns  for  the  avenge  and  the  variance  r.ignt 
be  worth  computing  roughly,  they  are  unlikely  to  make  enough 
of  a  contribution  to  bother  with. 

13  *  Numerical  Examples 

Let  us  consider  the  second  example,  and  suppose  that 
the  elements  to  be  used  are  5$  resistors,  iQjg  resistors  or  20% 
resistors  (i.e.  have  these  limits  on  deviations  from  their  nomi- 
nal  values)  and  let  us.  for  simplicity,  assume  that  their  dis¬ 
tributions  (within  these  limits),  are  uniform.  Then  the  necessary 
numerical  values  will  be  as  in  Table  3,  which  also  gives  the  for¬ 
mulas  for  the  cumulants  of  the  attenuation  in  terms  of  a 2  for  a 
nominal  attenuation  of  10.  These  formulas  will  converge  most 

slowly  when  all  resistors  are  20%  resistors.  In  this  situation, 
combining  teiros  of  the  same  order, 

ave  z  =  10.0000  +  0.1091  +  0.0029  +  ...  =  10.LU0  +  ... 
var  z  =  I.3389  +  0.0524  +  ...  =  I.3913  +  ... 

ske  z  =  0.0000  4-  0.2341  +  ...  =  0.2341  4- 

elo  z  -  -  1,0755  +  . . . 


mKmmmmmrnm 
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Table 


Numerical  quantities  for  uniformly  distributed  resistances  within 

specification 


Limits  Limits 


Resistor 

on  W 

on  w 

for  w 

y 

r 

G 

5# 

±.05 

±.5 

.0833 

0 

1.8 

0 

10$ 

±.10 

±  1 

•  3333 

0 

1.8 

0 

20$ 

±.20 

±  2 

1.3333 

0 

1.8 

0 

Formulas  for  cuniulants  of  attenuation  in  terms  of  oc  (through 


terms  in  c  ) 


ave  z  =  10.00000  +  0. 04091  {o^+oc?’)+0. 00074  (aA4+Op4) 


‘-.C'VI?  r  V/  +  . 


A  WC 


var  2  =  0. 16736 (a^  +  4aB2+oc2)+0. 00736( aA4+ac4) 


ske  z  =  -  0. 01642 (aA4+ac4)+0. 07381 (aA2+ aP2 )aR2+0. 01688 aA2aP2 


JA  '  “C  '  B 


+  ... 


elo  z  =  -  0. 03368(0/  +  l6aB4oc4)  +  ... 


Clearly  the  leading  terms  in  the  average,  variance  and 
elongation  need  little  adjustment.  The  first  term  in  the  skewness 
is  zero,  so  that  the  leading  correction  term  contributes  substan¬ 
tially.  There  is  no  reason  to  believe  that  the  neglected  terms  in 

1 

cr*  and  higher  will  affect  more  than  the  6th  significant  figure  in 
the  average  or  the  4th  in  the  variance.  If  this  be  so,  then  un¬ 
considered  terns  are  utterly  negligible. 


-  r  -'V  — «  V'-,'- ,v-:: ' . -  - 


^  v  ,1  #.-7  ^p 


111-12  - 


, .i  -  treatment  given  assumed  a  uniform  distribution 

of  resistance  within  toiler  .nc. a.  T’us  is  not  likely  to  be  the 

case  wif'  my  actual  simply  of  resistors.  In  practice  it  is  not 

unlikely  that  we  should  know  enough  about  the  actual  distribution 
to  make  a  somewhat  more  reasonable  assumption  but  our  knowledge 
will  always  be  far  from  perfect.  The  assumDtion  that  the  distri¬ 
butions  of  conductances  for  the  shunt  elements,  and  the  distribution 
ol  resistance  for  the  series  element,  are  all  uniform  within  tol¬ 
erances  is  manifestly  arbitrary.  But  it  is  not  a  bit  more  arbitrary 
than  the  previous  assumotion.  (We  might  exnect  the  logarithm  of 
resistance  or  ol  conductance  to  be  rlearly  uniform  under  certain  con¬ 
ditions.)  If  we  make  this  assumption,  then  we  have  to  deal  with  the 
numbers  and  lormulas  of  Table  4.  Here,  if  a  -  10  and  all  resistors 
are  20^  resistors: 


ave  z  =  10.0000  +  ...  =  10.0000  +  ... 
var  z  =  1.33*9  +  0.004*  +  ...  *  1.3437  +  ... 
ske  z  =  0.0000  +  0.2925  +  ...  -  0.2925  +  ... 
elo  z  =  -1,07547  +  . . . 
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Table  4 

Naiaei  leal  qu«nti  Litra  fur  uniform  cilai  ribuUun  within  tolerances 

for  shunt  conductances  and  series  resistances  erances 


Resistor 

Limits 
on  V 

Limits 

on  v 

c2 

for  v 

for  v 

y 

or  V 

P 

5Jf 

10% 

20% 

±0,05 
±0.1 
±0 ,2 

±.5 
±  1 
±  2 

0  OS) } 

0  3333 
1«333 

0 

0 

0 

1 

1.3 

1.3 

1.3 

0 
0 
0 

wh^Ta*-  Iff  Cumulants  of  attenuation  in  terms  of  a2  (of  v*s) 


ave 

var 

ske 

elo 


Z 

z 


*  10.00000  +  ... 

-  0.16736(orA2+4<TB2+crc2)  +  0.00136(<rA2+t-c2)  0 . 
“  °°7394  (aA2+<yc2)<TB^  .  0.01664  cfa*  *  ... 

■  -  0.03361(^+16 +  ... 


The  propagated  values  of  average  and  variance  are  some¬ 
what  different  from  the  results  obtained  under  the  previous  assump- 
tion.  These  differences  must  be  slight  from  a  practical  point  of 
view  since  they  derive  from  assumptions  which  are  practically  in¬ 
distinguishable.  Practical  accuracy  will  not  be  understated  if  we 
predict 


ave  z  *  10.0, 

var  z  -  1.35, 


ske  z  «  0,3  s 

elc  z  -  -l.l 


for  the  result  of  using  20 %  resistors  for  all  three  elements. 

Qualitatively,  the  indications  are  that  one  correction 
term  foi  t+.e  skewness  and  no  correction  terms  for  the  other  cumu¬ 
lants  (average ,  variance  and  elongation)  will  suffice  if  we  work 
with  suitably  transformed  individual  variables . 


I 

I 

$ 

| 


* 
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14.  Probabilities  of  deviations  beyond  tolerances 

Let  us  continue  the  previous  example,  supposing  that 
attenuations  between  8*5  and  12  0  are  acceptable,  while  those  out¬ 
side  these  limits  are  not  acceptable .  net  us  begin  by  calculating 
probabilities  by  the  normal  approximation.  For  all  resistors  20$, 
ave  z  -  10*12  and  var  z  =  1.34,  so  that  the  limits  lie  at 

standardized  deviates  of 


,  .  !.29  and  12. °  ,10^  -  1.72 

y/TTW 


respectively.  Reference  to  a  table  of  the  standard  normal  distrl 
bution  indicates  probabilities  of  9.85$  for  too  low  attenuation, 
and  4.27$  for  too  high  attenuation. 

If  we  wish  to  make  better  approximations  we  find  first 


that 


?1  = 


0.3 

(l.35)3/^ 


0.19 


72  = 


-1.1 

(1.35)2 


0,60 


and,  consequently,  =  0,036,  @2  =  2.40.  Reference  to  the 
Merrington-Pearson  table  (specifically  table  42  in  Pearson  and 
Hartley)  for  these  values  of  3^  and  3g  (and  positive  skewness) 
leads  to  the  following 


, h' •>  ..  H  .  ■,  ■ 

glsgg  :  ■  7m  ■ 


<*w&-  \  .  *  .  i  *. 
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Lower  Deviate 

Normal  .036,2,40  Diff. 

Upper  Deviate 
Normal  .017,  2.78 

Diff , 

5# 

-1.645 

-1,57  +.075 

4-1.645 

1.72 

4- .  07^ 

2.5* 

-1.960 

-1.76  4-.20q 

4-1.960 

1.99 

+  .03n 

1# 

-2.325 

-1.93  4-.395 

4-2.325 

2.26 

-.065 

We  are  thus  led  to  modify  our  standardized  deviates  as  follows 

-1.29  +  0.01  =  -1.28 
+1.72  4-  0,06  =  I.78 


These  modified  deviates,  when  looked  up  in  a  unit  normal  table 
correspond  to  10. 03#  below  and  3.75#  above,  the  changes  in  the 
two  tails  tending  to  compensate  one  another.  The  total  #  out¬ 
side  tolerance  decreases  by  O.52#  -  0.180  =  O.34#  out  of  about 
14#.  Clearly  this  change  has  little  practical  significance. 

Thus,  as  illustrated  above,  the  methods  of  using  higher 
cumulants  are  often  of  most  value  when  they  demonstrate  the  un¬ 
importance  of  the  correction. 


I 
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TRANSFORMATION  if  TERMS  OF  RESPONSE 


15.  General  formulae 

We  have  now  to  consider  the  possibilities  of,  and  neces¬ 
sary  machinery  lor,  changes  In  the  terms  In  which  the  response  is 
measured.  We  shell  clearly  wish  to  keep  partial  effects  linear. 
This  requires  us  to  change  the  terms  of  the  Individual  variables 
at  the  same  time  that  we  change  the  terms  of  response. 

If  z  =  g(v1,v2,  .  .  . ,  vk)  was  such  that  g&a,  gaa&,  etc. 
all  vanished,  we  want  to  apply  y  =  <p(z)  and  compensate  for  this  by 
also  applying  vfl  =  va(ua),  so  that  In  the  resulting 


y  f (u^,u2>  =  •  .,  u^. ) 


we  shall  have  ail  ffla,  f&aa,  etc.  vanish.  If  we  maintain  the 
verbal  definition  of  the  —  the  standard  deviation  of  the  partial 
effect  of  individual  variable  a  -  rather  than  maintaining  its  nu¬ 

merical  values,  so  that  now 


T 


a 


2 


var  u  , 
a 


then  the  propagation  formulas  Involve  the  xa  and  the  reduced  higher 
(cross)  derivative  values  of  f,  defined  by 


s 


'ab 


ab 


f  f, 
a  b 


s 


aabbb 


faabbb 

"^7T~ 
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Sj»d  the  lilcc.  W«?  c  all  .  f 'J’f'.ra  *•  •->  1  w  ^  v  ,  ,,,  , .,  ,  , 

'  Iorr  lately,  express  the  s*s  rather  simply 

In  terms  of  the  *  *s  and  the  derivative  values  of  *, 


v  -  c)  , 


(~ §-> 

dz  '  ( aye  z ) 


The  resulting  formulas,  derived  in  part  VIII, 


are 


Sab  ^  +  *" 


abb  ~  ‘abb  +  ‘ab  +  1>"  '  -  <P>", 


!  _  ,,  M 


5abbb  ~  tabbb  +  (2<P"’  -  3 <pV)  tab  +  cpiv-  V'V  +  3(pyy.j 


~  t---’-  +  <P"  (t  .  +  2  t  2 


'aatb  aabb 


aab  '  *“  vab  +  t  .  .  ) 

abb' 


+  (V"  -  3<p"cp")  tab  +  Cp:LV  -  29"V  +  9yy 


and,  indeed. 


saaabb  “  ‘aaabb  +  6t?"  ‘ab^'aab  +  -  <P'V')taab 

+  (29"'  -  39’V1)  tabfe  +  6  ( qp " '  -  9'V')tab2 
+  (6fv  -  l6cp"iq,"  +  9<p "jp 'V')t  b 

+  9V  .  W'  -  <p"V"  +  Jcp"  y9 "  .  39'VV'ip " 

•  i*vf- .  .  ••  -  -  ** ;  ^  '■■  'ywy  pgg?»ft»<yv. 

'  ••  '  *'  ■  ,Z-  '■££  »  't-‘  ,,  ”C  '  .  l3  y*, 

*.  >'r.i '■’•.■•  ••'  v  **  :'£-:ry,.  ,  '•■  •  ._>? '  fc  Jr  '  'j& 

■  V  •  r  •"*  •  *  ’  T;-  "f-  •  4 

' ,  -4-  ■  ■  : 
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* ^ »  Sped  a^lgaf  ions  for  powe r  transformations 

We  are  usually  not  going  to  use  a  more  complicated 
transformation  than  y  -  A(z+c)p  (or  one  of  the  limiting  forms  of 
this  family ,  which  include  logarithmic  and  exponential  forms). 

The  higher  derivative  values,  qp  can  be  thrown  back  onto 

<p',  q>"  and  <p"'  for  any  such  power  transformation  (see  Part  VIII 
^  or  details)  and  If  we  fix  q?  *  =  1  and  use  these  relations,  we  may 
find  fom  alas  por  directly  changing  the  actual  coefficients  of 
Table  2.  These  i irmulas  are  given  in  Table  y.  The  abbreviations 

L  =  <p"»  -  <pV' 

and 

E  =  cpn(j, 11 


contribute  substantially  to  compactness. 


Effect  on  coefficients  of  a  power  transformation 
(with  <p  *  =1)  of  response  and  related  transforma¬ 
tions  of  the  Individual  variables. 

Abbreviations:  D  *  c?,!t  -  <p'V'*  E  «  <p Hcp ” 


Vabb  =  taatb  +  ^"^aab  4  2tab  +  tabb^ 
+  (4l>fE)t  .  +  -At  D(2D+E)  , 


saaabb  fcaaabb  +  6<?"  tab  taab  +  6D  fcaab  +  (2D’E)  ti 


+  6D  tab?  +  r  (12I)2  +  2DE  •  E?)  lab 
+  |  (6lP  +  2D2E  -  DEC)  , 


8aab  +  sab  +  sabb  =  ^ab  +  fcab  +  *abb  +  V  *ab  +  (2D+E)  - 


"Sf  aaabb  +  3ab  aaab  +  3  saaab 


1  1 
2  ^aabb  +  ^ab  ^aab  +  3  ^aaab 

+  *"  ^ab  +  2tab2  +  W  +  5  (22»+7E)  t 
+  5^7,  (10D2  +  7DE) 
fcab  +  f" 


3saab  +  a  sabb 


3taab  +  I  tabb  +  I 


IV-5  * 


One  substantial  advantage  of  Table  5  Is  that  it  allows  us 
to  calculate  the  effect  of  a  power  transformation  with  chosen 
v'\end  <p,,f  on  a  particular  set  of  propagation  formulas  without  any 
necessity  of  calculating  the  powei  transformation  explicitly.  This 
often  suffices,  by  showing  the  most  that  can  be  done. 

Actually,  as  13  shown  by  one  of  the  examples  (see 
Section  37) >  the  apparent  form  of  the  power  transformation  is 
sometimes  shifted  very  much  by  small  changes  in  9"  and  9"'*  The 
nature  of  y  =  <p(z)  near  z  =  z^  will  not  change  much,  but  the 
analytic  form  may  change  considerably  —  as  from 

y  =  (400) ^  (z  +  -y-)1/3 

to 

y  *  -  2l5.7(z+i4.42)~-2757 

which,  apart  from  an  additive  constant,  are  quite  similar  near 
z  =  0.  The  c  in  (z  +  c)p  imparts  more  flexibility  for  a  fixed  p 
than  one  naturally  supposes. 

If  one  wishes  to  restrict  attention  to  the  simple  power 
transformation  y  =  kzP ,  then  only  9"  may  be  chosen,  and 

9";  =  9  "9"  -  <p"  ,  D  =  -  cp'/z  .  E  =  9  "9" 

o  /  o 

are  an  inevitable  consequence. 


# 
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27.  Made. -i  of  use 

We  may  use  this  freedom  of  transformation  of  terms  of 
response  in  various  ways,  the  most  promising  of  which  seem  to  be 
either  further  simplification  of  formulae  or  reduction  in  skewness 
of  response.  More  detailed  inquiry  in  Part  VIII  leads  us  to  the 
views  (l)  that  the  latter  will  only  infrequently  be  practical, 
since  its  effective  accomplishment  depends  on  detailed  quantitative 
knowledge  of  the  and  xa  and  (ii)  that  effort  on  simplification 
cf  formula  is  best  concentrated  on  the  annulling,  or  reduction  in 
size,  of  the  coefficient  of  r  2  2  In  the  variance. 

a.  D 

Since  this  coefficient  is,  for  any  particular  (a,b) 


3aab  +  sab  +  sabb  “  fcaab  +  tab2  +  tabb  +  V'  tab  +  V  "  <pV 


and  since  we  have  both  qp"  and  cp'1  *  at  our  disposal,  we  can  exactly 
annul  two  different  coefficients  by  a  proper  choice  of  <p.  (How 
effective  such  a  transformation  may  be  will  depend  very  much  on 
the  particular  example.) 

The  examples  considered  later  on  indicate  that  the 

apparently  simpler,  and  certainly  less  justified,  approach  of 

trying  to  annul  some  of  the  sab  alone  is  not  likely  to  be  rewarding 

If  we  are  to  gain  in  simplicity  of  formula  by  transforming  the 

response,  we  must  consider  third  derivatives  as  well  as  second 

derivatives  so  as  to  obtain  and  use  t  and  t 

aab  abb  * 


If  we  consider  transforming  the  response  In  the  first 
example ,  we  find  that  we  can  eliminate  the  main  correction  term 
the  formula  for  the  variance.  Instead  of 

var  z  £  Z  +  0.01  Z  t21_12  t212 


we  obtain 


var  z  Z 


2  Ta  +  °-00°5  s(>2i-lT2i-l+7 


21t21 


r2i-l 


21 


If  we  had  needed  to  reduce  this  term,  the  gain  would  be  striking, 
even  though  the  simplest  transformation  uses  the  cube  of  the 

difference  between  the  time  delay  and  a  constant  larger  than  any 
likely  delay. 

Transformation  of  the  response  In  the  second  example 
again  reduces  the  first  correction  term  In  the  variance  by  a 

factor  of  50  or  so,  but  again  this  reduction  was  not  badly 

needed.  The  power  transformation  which  does  this  Is  proportional 
to 


interestingly  enough. 
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V  DETAILS  OF  PROPAGATION 


Cocumulants  -  definitions  and  properties. 

The  covariance,  defined  by 

cov  (x,y)  »  ave  xy  -  ave  x  ave  y  =  ave  (x-x)(y-y) 

and  satisfying 


(i) 

var  x+y  =  var  x  +  2  cov  (x,y)  +  var 

y 

(11) 

var  x  =  cov  (x,x) 

:| 

(ill) 

cov  (x, y)  =  cov  (y,x) 

3 

(lv) 

if  y  is  (statistically)  independent 

of  x,  then 

cov  (x,y)  =  0 

(v) 

var  2xa  =  2a2b  cov  (x^) 

M, 

'll 

(vl) 

cov  (2xa,  2yb)  =  2^  cov  (xa,yb) 

(vii) 

if  v  is  independent  of  x ,  y  together,  then 

cov  (x,  yv)  =  cov  (x,y)  ave  v 

.  K 

"w 

has  stood  beside  the  variance  for  a  long  time.  Similar  partners 
for  the  higher  cumulants  do  not  seem  to  have  been  provided,  though 
they  will  prove  very  convenient  to  us. 

We  may  define  the  co skewness,  a  function  of  three 
arguments,  by 

cok  (x,y,z)  =  ave  (x-x) (y-y) (z-z) 


t 


(we  use  "cok”  because  "cos "  Is  preempted)  and  can  show  directly 
that 

(l)  ske  x+y  -  ske  x  +  3  cok  (x,x,y)  +  3  cok  (x,y,y)  +  ske  y  , 
(ii)  ske  x  =  cok  (x,x,x)  , 

(ill)  cok  (x,y,z)  -  cok  (x,z,y)  -  .  .  .  .  =  cok  (z,y,x)  , 

(iv)  if  z  is  (statistically)  independent  of  x,y,  together, 
then  cok  (x,y,z)  *  0 

(v)  ske  2xa  =  2a£fa2c  cok  (x&,xb,xc) 

(vi)  ook  (Zxa,  Zyb,  Z2c)  =  2aZbS0  ook 
( vi I )  If  v  Is  independent  of  x,y,z  together,  then 

cok  (x,y,zv)  =  cok  (x,y, z)  ave  v 

in  entire  analogy  with  the  properties  of  the  covariance. 

Moreover,  we  can  define  the  coelongatlon,  a  function  of 
4  arguments,  by 

coe  (x,y,z,w)  =  ave  (x-x) (y-y) (z-z) (w-w) 

-  cov  (x,y)  cov  (z , w) 

-  cov  (x,z)  cov  (y,w) 


-  cov  (x, w)  cov  (y,z) 


-  v-3  - 


and  can  show  direct:,  that 


(i)  elo  x+y  -  elo  x  +  <♦  coe  (x,x,x,y)  +  6  coe  (x,x,y,y) 


+  4  coe  (x  .  V  .  v  ir  ^  4.  al 

%  '  W  >  */  / 


(ii)  elo  x  =  coe  (x,x^x,x) 

(ill)  coe  (x,y,z,w)  =  coe  (x,y,w,z)  =  .  .  .  =  coe  (w,z,y,x) 

(lv)  if  w  is  (statistically)  independent  of  x,y,z  together, 
then  coe  (xj^z^w)  =  0 

(v)  elo  2xa  =  £a2b202d  coe  (xa,  x^,  xc<  X(J) 

(vl)  coe  (2;:a,  2yb>  S«e,  2W(J)  -  ZaWd  coe  (xa,  yfa,  Z(j,  „„) 

(vli)  if  v  is  independent  of  x,y,z,  w  together,  then 
coe  (x,y,z,wv)  -  coe  (x,y,z,w)  ave  v 

again  in  entire  analogy  with  “le  properties  of  the  covariance. 

20 •  -Explicit  formulas  in  two  forma 

By  using  x-x,  y-y  etc.,  we  have  kept  the  defining  equa- 
tlons  for  the  cocumulants  moderately  short.  There  will  at  times 
be  something  to  be  gained  by  writing  these  expressions  out  in 
detail.  They  appear  as  follows; 


-  V-4  - 


c  jk  v  y* z )  *  ave  xyz  -  (ave  x)(ave  yz)  -  (ave  y)(ave  xz) 
*(ave  z)(ave  xy)  -f  2(ave  x}  (ave  y)(ave  2) 

*  ave  xyz  -  (ave  x)(ave  y)(ave  z) 

-(ave  x)  cov  (y,z)  -  (  ave  y)  cov  (x,z)  - 

(ave  z)  cov  (xsy) 

coe  (x,y,z,w)  -  ave  xyzw  -  (ave  x)(ave  yzw) 

-  (ave  y ) (ave  xzw)  -  (ave  z)(ave  xyw)  -  (ave  w)(ave  xyz) 
(ave  xy)(ave  zw)  -  (ave  xz)(ave  yw)  -  (ave  xw) (ave  yz) 

•+  2(ave  x)(ave  y)(ave  zw)  +  2(ave  x)(ave  z)(ave  yw) 

+  2 (ave  x) (ave  w)(ave  yz)  +  2(ave  y)(ave  z)(ave  xw) 

+  2 (ave  y ) (ave  w) (ave  xz)  +  2(ave  z)(ave  w) (ave  xy) 

-  6(ave  x) (ave  y)(ave  z)(ave  w) 

-  ave  xyzw  -  (ave  x) (ave  y)(ave  z)(ave  w) 

(ave  x)  cok  (y,z,w)  -  (ave  y)  cok  (x,z,w ) 

(ave  z)  cok  (x,y,w)  -  (ave  w)  cok  (x,y,z) 

cov  (x,y)  cov  ( z j w)  -  cov  (x,z)  cov  (y,w)  ( 

-  cov  (x,w)  cov  (y,z) 

(ave  x)(ave  y)  cov  (z,w)  -  (ave  x)(ave  z)  cov  (y,w) 

-  (ave  x) (ave  w)  cov  (y,z)  ~  (ave  y)(ave  z)  cov  (x,w) 

-  (ave  y) (ave  w)  cov  (x,z)  -  )ave  z)(ave  w)  cov  (*,y) 

The  apparent  lack  of  symmetry  In  the  basic  definitions 

can  be  removed  by  rewriting  them  In  the  forms; 


V-5  - 


ave  x  *  ave  x  , 

ave  xy  -  (ave  x)(ave  y)  +  cov  (x,y), 
ave  xyz  =  (ave  x) (ave  y)(ave  z) 

+  (ave  x)  cov  (y,z)  +  (ave  y)  cov  (x,z) 

+  (ave  z)  cov  (x,y)  +  cok  (x,y;z), 
ave  xyzw  =  (ave  x)(ave  y)(ave  z)(ave  w) 

+  (ave  x)(ave  y)  cov  (z,w)  +  ave  x)(ave  z)  cov  (y^w) 

+  (ave  x) (ave  w)  cov  (y,z)  +  (ave  y)(ave  z)  cov  (x,w) 

+  (ave  y ) (ave  w)  cov  (x,z)  +  (ave  z)(ave  w)  cov  (x,y) 

+  cov  (x,y)  cov  (z,w)  +  cov  (x,z)  cov  (y,w) 

4'  cov  (x,w)  cov  (y,  z) 

+  ave  x  cok  (yjzjw)  +  ave  y  cok  (x,z,w) 

+  ave  z  cok  (x,y,w)  +  ave  w  cok  (x,y,z) 

+  coe  (x,  y,  z,  w) . 


Extension  to  higher  order  cumulants  is  now  obvious. 

These  formulas  appear  messy,  although  they  are  really 
quite  simple.  A  different  notation,  which  Is  not  as  convenient  for 
working  with  covariances,  coskewnesses  and  coelongations  as  tools, 
shows  this  more  clearly  We  use  it  here  for  Immediate  expository 
use  only.  Put  ave  x  =  x,  cov  (x,y)  =^77  ook  (x,y>z)  .  £7^ 
coe  (x,y,z,w)  =  x,y,z,w,  then 


-  V-6  - 


x,y  *  xy  -  x  y 


x,y,z  *  xyz  -  x  yz  »y3c2-?5cy+?xyi 
*  xyz  -  Z*  a  be  +  2x  y  z 

x,y,"z,w  *  xyzw  -  2*  a  Ec3  +  22*  ab  cd  -  6x  y  z  w 


xy  =  x  y  +  x,y 


xyz  =  xy  z  +xy,z+yx,z+  z  x,y  +  x,y,  z 


-  x  y  z  +  2  a  b,c  +  x,y,z 


*  _ 


xyzw  -  x  y  z  w  +  2  a  b  c,d  +  2"  a,b  c,d  +  2*  a  b,c,d  +  x,y,z,w 


where  2  has  the  same  interpretation  as  elsewhere  in  the  memorandum, 
and  a,b,c,d  is  some  permutation  of  x,y, z,w.  The  formulas  are  much 
more  perspicuous  in  this  notation.  (if  the  notation  were  to  be 
extended,  var  x  =  x,  ske  x  =--  x,  elo  x  =  x  might  be  used. ) 

21.  Independent  centered  monomials 

We  next  consider  the  cocumulants  of  expressions  of  the 

f°rm  waV  where  wa  and  wb  are  (statistically)  independent  with 
moments 


0,  or 


3 

7  0, 

'a  a 


r-°-  Ga°a5  and  °' 


a”a’ 


n°b ' 


Vb 


a 


V 


respectively.  Using  the  available  relations,  these  can  be  easily 
reduced.  Thus,  for  example. 


. . 


7-7 


O 

I 


3  2 

ave  wa'wb 

-  ave  wr wK  ave  w, 
a  0  e 

3 

ave  wv  ave 

C2 

wb2  -  (ave  wa)(. 

„  2 

>„a,  a, 

7  a  a  b 

*  *  ♦ 


since  ave  v/a  »  0.  Similarly,  dropping  at  once  terms  which 
clearly  vanish, 

cok  (^a,wawb,  wa2wb2)  =  ave 

2  2 

-  ave  wa  wb  cov  (*a>wawb^ 

*  VaW  -  (ffaV,!0) 
■  VaW 

By  similar  calculations  we  reach  the  results  shown  in  Tables  6 
and  7,  as  well  as 


coe 

(wa 

>wa 

>wa 

-V  - 

<V 

3)aa 

4 

coe 

(wa 

,wa 

’wa 

•V  * 

0 

coe 

(wa' 

,wa 

’wa 

,wa2)  -  (Ga 

-^a 

coe 

iwa, 

,wa 

•waV 

*  0 

coe 

(wa, 

,wa 

>wa 

•«b2)  ' 

•  0 

v-o 


TABLE  6 

Table  of  coefficients  by  which  cr  *",’m  a,  °+n  must  be 

cl  D 


multiplied  to  obtain  the  covariance  of  w  1  w.  u'  with  w  *“  w  “  when 


a 


m  ,  n 
v/, 

a  b 


wa  and  w^  are  (statistically)  independent,  and  each  averages  zero, 


w_ 


a 


in 


W 


a 


W  W. 

a  b 


v/. 


w 


a 


w. 


w 


a 


w  w, 
a  b 


w. 


w_ 


a 


a.  b 


W  W. 

a  b 


w, 


w. 


3 


3 

w0Jw. 
a  b 


2..  2 
b 
3 


a  b 


w  wb 
a 


wx 


0 

7, 

C 

0 

0 


0 

1 

0 

G, 

< 

0 

y 

*b 

o 


a 


0 

1 

0 

0 

?b 

0 

1 

0 

rb 

0 

?a 

?b 

0 

Gv. 


0 


V1 


0 

0 


G  -y 
a  7  a 

0 


0 


V1 


Vb 


0 


0 

0 

0 

1 

0 

0 

ya 

*b 

0 

0 

Fa 


Vb 


0 


0 


0 


V1 


0 


0 


Gb-?b 


0 


^a^b 


V1 


0 


0 


G  -r 
a  7  a 


0 


Note.  With  appropriate  subscripts  var  w  =  2 ,  ave  w3  =  7a3, 

ave  w^  =  To^f  ave  w5  =  Go5. 


V-9  - 


TABLE  7 


Coefficients  of  a  1+m+p<T.  ^+n+G 


in 


cok  (wa1eb  i,  w.VU.**,*) 


_  r>  -i_  J  rt  vv-i 

a  b  a  b 


When  wa  and  wb  are  (statistically)  independent  and 

each  averages  zero. 


1-2  2  J  %  wa2  wawb  "h®  »- 3 


a 


w„  w 
a  a 


ra  0  ra-i  0 


wa  wb 


ooo 


0 

0 


G  -  7 
a  'a 


0 


Wb  wb 


0^0 


VI  0 


wa  wa  V1  0  V2?a  0 


0 


wo  w  w. 
a  a  b 


0 


1  0  7-1 
'a 


wa  wb 


w,  w 
b  a 


0  0  0 
0  0  0 


?b 

0 

0 


Wb  wawb 


0 


0 


?a  V1 

0  V2?b 

— Wlth  appropriate  subscripts 


wb  wb 


0  rb-i  0 


var  w  =  c  , 

ave  vr^=  yo^ 

4  4 

ave  w  =  ro  , 


ave  Ga*. 


2 

wa  wb 


0 


a 


2  3 

wawb 


a 


?b 

0 


0 


0 


G,  -v, 
b  '  b 


-  V- 


1  A 

A  v 


22‘  I!5Por>  styles  In  Independent  guanti 
Suppose  that 


f*  f  s»  £* 

V,*,  w  5? 


2  «  h(w1,w2/...,wk) 

where  ave  wfi  =  0,  and  the  wfl  are  Independent.  Then  the  usual 
multiple  Taylor  Series  for  «  can  be  written  m  terns  of  the  derive- 
tlve  values  at  (0,0 . 0),  which  Is  the  average  point,  such  as 


aab 


dwa  dwb 


^ C *  ^2 * 1  •  •  ) 


0 ,  ,  .  .  j  0 ) 


In  carrying  out  this  development  we  will  make  full  use  of  £*,  a3 

used  above  and  discussed  In  Section  38.  When  we  do,  the  Taylor 
series  becomes 

z  -  h(0,0,...,0) 

+  Sh  w 
a  a 

+  2  Shaawa2  +  2\bwawb 

+  5  Zhaaawa3  +  §  ^\ab«a\  +  jV.»A« 


abc  a  bwc 


+  ®  Zhaaaawa4  +  §  2\aab«a\  +  5  **W«V 

+  2  2  haabcwa  wbwc  +  2  habcdwawbwcwd 
+ 

'  •  *  * 

Where  the  symmetry  of  h  and  h 

aabb  a  anc^  b  is  to  be  recognized 
in  interpreting  S  ,  and  where  all  the  coefficients  can  be  easily 

checked  by  direct  differentiation- since  each  monomial  appears  once 
and  only  once. 


■  -  |F 


-  »«( V  .  • 


7-11 


From  this  expansion,  slightly  extended,  we  Immediately 


deduce 


f 

w 


ave  z  =■  h(0,0,...,o)  +  |  Zhaa<>a*  *  5  ^*aaa  7aoa3 

1 
7 


,  -  4  1  * 

+  A-i,  Zh _ r„o„  +  4  2  0  "c  "  + 


1 

24 


aaaa  a  a 


2  2 

‘aabb  ua  ab 


TZT-  £haaaaa°aca’ 


+  T5  Z\aabb  >aV  V  +  • 


3.  2 


If  we  think  of  writing  out  the  variance  of  z  In  terms  of 
a  double  sum  of  covariances  of  its  terms,  and  use  Table  6,  we  see 
that  we  have 

2  2 

var  z  =  2h  a 
a  a 


+  Zh  h  y  a  ' 
a  aa'a  a 


+  5  aahaaara0a4  +  *  “a^V1)®,4 
+  SXhabb+hab2+haabV  °aV 
+  h  ^VaaaaVa5  +  \  aaahaaa<  W aa5 


+  S  (2hahaabb  +  2haahabb  +  habhaab  +  IhaaabV  Va^ 

+  terms  of  order  >  a 

where  the  symmetry  of  tht#whole  coefficient  of  c&2cb2  is  to  be  rec¬ 
ognized  in  Interpreting  2 — a  not  unmanageable  expression. 

Writing  ske  z  as  a  triple  sum  of  coskews  of  the  terms  of 
z,  and  using  Table  7,  we  find 


-  V-15  - 


ske  z  m  Zh  '‘ye* 

H-  <£ 


2 _  2 


+  2  Zha  hzz^a'1}  aa~  +  6  1  Wab  0a  a] 

+  5  ffla2haaa(oa*’'a)oa5  +  I  ’Vaa^V’VV 

+  ^Vab2^.-1)  ».  V 


a  b 


+  2*  (^-ha2habb  +  3haahabhb  +  3hahaabhb)^fla3ab2 


+  terms  of  order  >  a  . 


The  same  approach  to  elo  z  leads  to 


elo  z  -  ZhaVa-3)  ®a4 

fh 

+  terms  of  order  J>  a  . 

23 •  The  generalized  propagation  formulas 

The  formulas  in  the  last  section  were  derived  on  the  , 
apparently  equally  important,  assumptions  that  the  w  were  inde- 

St 

pendent  and  had  average  zero.  Suppose  them  independent,  but  their 
averages  not  to  be  zero,  and  put 


°a  =  wa  *  ave  wa  ' 

then  the  qa*3  will  satisfy  both  assumptions,  and  the  formulas  will 
apDly  in  terms  of  the  q&,  provided  we  expand  around  q^  =0,  i.e. 
wa  =  ave  w  .  Thus  if 


;■ 

I 

1 


M 


v-13 


V 


ca  *  ave  *"  aV€  (w  -ave  w  , 

3*  Jfc 

3 


^aaa  c  ave  =  ave  (wa-ave  wo) 


*a 


a 


4  jf 

F_a  *  ave  q  *  ave  (w  -ave  w 

ri  '3  «a  * 


a  a 


C  jr 

Gaaa  =  ave  =  ave  (wfl-ave  wo)- 


and  if 


ab 


£w  dw,  ^(w2'w2>  *‘*'wk:) 

a  d  K 


(ave  w.p  ave  Wp,  . . . ,  ave  xv,  ) 


and  so  on,  then  the  formulas  for  var  z,  ske  z,  and  elo  z  will  hold 

without  change.  A  small  change  will  be  needed  in  the  formula  for 
ave  z,  and  it  will  become 


ave  z 


h(ave  w^,  ave  w?,  . . . ,  ave  w^J 


+  i  Zh  a  2 
2  aa  a 


+  ■T'  Zh  y  o  ^ 
o  aaa  7a 


aaa  'a  a 
1  4 


2  5  Zhaaaa  raGa  +  If  2  haabb  aa2a'  2 


a  b 


+  Tfc  Zhaaaaa  Va5  +  h  Z\aabb  V 
+  terms  of  order  >  ft . 


1 


•I 


>. 


and  so  on. 

We  wish  to  Illustrate  the  effect  of  making,  or  omitting 
various  transformations.  But  we  would  avoid  needless  complexity. 
So  we  shall  treat  the  case  where  each  section  has  the  same  design 

value  in  considerable  detail,  and  then  treat  the  general  case  by 
the  quicker  methods. 

If  all  sections  are  to  be  alike,  we  may  select  units  of 
Inductance  and  capacitance  such  tha.  she  design  value  of  each 
element  Is  10.  This  means  that  ca  =  1  for  a  component  with  a 
standard  deviation  of  10*  of  its  nominal  value.  Hence  values  of 


-  VI -2  - 


Ca  as  large  as  2  are  possible  but  not  unlikely,  while  values 

between  1  and  9.1,  or  possible 

'  *ghtly  2mail  —  may  be  regarded  as 

likely.  This  means  that  we  can  draw  relatively  fair  Impressions 

from  numerically  stated  formulas  Involving  different  powers  of  „  . 
With  this  choice,  we  find  a 


t 

4* 


:?S 


ha  "  0-5, 
haaa  =  00375, 


Jaa  ~  -  0.025, 


r‘aaaa  ~  -  0*00094 


h 


aaaaa  O.OOQ33  , 


ab 


’0.025, 

:o, 

•0.00125, 


aap 


aabb 


h 


aaab 


(0, 

[o,,oooo6+, 

> 

► 

► 

I  - 

0, 

■0.00019, 

(C, 

’-0.00001, 


h 


aaabb 


0, 


if  (a,b)  -  (21-1,21)  or  (21,21-1) 
otherwise 

if  (a,b)  =  (21-1,21)  or  (21,  21-1) 


otherwise 

if  (a,b)  =  (21-1,2!)  or  (21, 


21-1) 


otherwise , 

(a,o)  -  (21-1,21)  or  (21,  21-1) 
otherwise 

if  (a,b)  =  (2i-i,  21)  or  (2i,  2i-i) 


otherwise 


VI -3  - 


hence 


ha  -  0,25,  hahfia  .  -0.0125,  hahaaa  -  0. 00187+ 
.+ 


h--  "  °-ooo6i  »  h-Ji - -  -0. 00047' ,  haahaaa  -  -.00009+ 


aa 


a  aaaa 


ha  "  °-li:5‘  hahaa  "  -0.00625,  hfh  =  0.00094',  h  h  2 


ha  -  0,0625,  h^haa  *  0.00312+ 


a  aa 


O,0003l7 


and,  provided  (a,b)  =  (21-1,  21)  or  (21,  21-1)  aiso 

^a^abb  +  !lab  +  ^aabhb  =  '  °.00062+ 

1  1 

f  hahaabb  +  2  haahabb  +  habhaab  +  5  'Wb^b  =  0  -00003 

6hahahab  "  °-°375 

3hahab  "  0.00094" 

3  2 

2  hahabb  +  3haahabhb  +  3hahaabhb  =  -O.OO236, 

theae  expresalona  vanlehing  for  all  other  palra.  We  can  now  write 
down  the  propagation  forraulaa,  finding 


ave  z  -  10 j  -  0.0125  2a 


a 


+  0.00062  Zy  a3  -  0.00004  ir  a^ 
a  a  a  a 


+  0.00002  2o sl,fosf 

+  0.00000  26  0?  -  0.00000  2 (>„,  ,0  3-  2  .  v  „  2  3 

8  8  l72i-l°21-2a2i  +  721 °21-la2l) 


0  c  0 


71-4 


t 


var  2  «  0.25  Zo 


a 


-0.02250  Zy c|  +  0,00062  zr  a4 

a  a  a  a 


a  a 


+•0.00016  2(ra-l)0g  -  0  00062  £o31Jc22 


-0.00004  SO  <rj  +  0.00002  2(0  -v  )05 
a  a  v  a  'a'  a 

+0.00003  2(y2i-la2i-la2i  +  r2ia2i-la21 ) 


ske  z  =  0.125  2  7  3 

a 


-  O.OO937  2(ra-l)aa  +  0.0375  Sa21_^o2^ 


+  0.00047  2(0a-7a)<^  +  0.00023  2(0-2?  Jo5 

a  a  a 

+  0.00094  +;(,81-l)a81.f<%3) 

-  0.00236  2(?21.1c21_3o22  +  ?21o21_j023 


elo  z  .  0.0625  s(ra-3)aa4  +  0.00625  2(^-4?^  +  .  .  . 


where  the  index  a  rune  from  a  to  2j  and  the  index  i  runs  from  1 

to  J.  We  can  simplify  these  somewhat,  with  the  results  already 
quoted  in  Section  10. 


**£*>£•«*** a*HW  •waseewMStt 


;  if  $$?hj 


1  *  .  .-1- 


’ 


if  l~ '  ^  * 

»-i 


f  -  . ..  -4-  a.V.^'.is%S-3!5--V. 


' ~  *  „  ■  *  '  ft-  :  ••  ■%■  JL".  .V'C:\ 

,  '  ,  ' 

, 


&  -r*  *.  if  ft  Jtfle.  s  «. 


25 .  T ranaforml ng  indiv 1 du ale 


If  now  we  wish  to  Introduce  v*s  whose  partial  effects  will 


he  linear,  we  have  only  to  look  at 


h(L1#  c2,  ,  .  *  >  L2 j_ i -»  C2 j )  ~  ^L1C2  +  •  +v/L2j-ll 


to  see  that  we  should  take 


-1  '  A  ^21-1' 


v2i  A  ^C21  ’ 


and  if  we  choose  A  so  that  the  first  derivative  value  is  unity. 


this  means 


.-1  =  2v/1°  V'/L2i-1 


v01  ~  2  yio  v/Cn<  * 


^L2i-1  C2i  "^O^i-l  V« 


whence 


Z  =  S(vi,v2,...,  v2  ) 


40  <V1V2  +  v3v4  +  •  •  •  +  V2j-1  v2j) 


’  *  wmm 


VI-6  - 


and 

S2i~l  =  £o  V2i  *  2 

g2I  =  roV2i-2  =  i  ' 


m 


fe 

ifjf  1  .  :  ' 

S2i-l,2i  ~io 

I 

i  a11  other  derivative  values  vanishing,  (we  note  for  later  use  that 

1 

^21-1, 21  =  0b1  all  other  tfs  vanishijig.) 

The  coefficients  which  appear  in  the  propagation  formulas 
when  (a,b)  =  (21-1,21)  or  (21,  21-1)  are,  omitting  those  which 
vanish, 

sab2  =  IWo  "  °’00O62+  , 

6sagbsab  =  TTO  =  8%  "  0-°375  , 

3®asab2  =  3KTO  “  °- 00094  "  > 


and  the  resulting  propagation  formulas  are 
ave  z  =  10J  +  no  other  terms 
var  z  *  0,25  2a2 

Si 


+  °'°375 

+  0.00094  S[(721.1-l)a21.;1+(r21-l)a21]021_2a2^ 

+  exactly  two  sets  of  terms  of  order  , 

elo  z  =  0,0625  £(r  -3)0^ 

9.  a 

+  a  finite  number  of  terms  of  orders  0 ^  to  . 

26-  The  general  case 

If  we  return  to  the  general  ease  of  a  delay  line  made  of 

possibly  unequal  sections,  we  have  (with  a  change  in  scale  for  the 
va  for  convenience). 


VI -8 


where  t±  is  the  tits  delay  contributed  by  the  1th  section,  all 
other  t  *s  vanishing.  The  resulting  formulas  arp 


t21-1  ~  (C2i )  var 


T2i  =  ^L2i-1  ^ 


ave  2  =  2  t1  , 


var  z  =  2  x 


T  2  2 
t21-1  t21 


ske  z  =  2  >  t  ^ 
'a  a 


6  m-  2  2 

4  Toj  T  T r\j 


i  21-1  t21 


l.^^ai-rn^^H-ty^-UTp,  ] t21  _ a 2 Tp  2 


+  terms  of  order  > 


T  ‘  '  '  ~~~  '  ‘''T’llr*  ’  k  ';V 

-  si 


>-l?fcrS:3C'‘,’5:  £ 
:■  *  ■ 


ESSa 


VI-9 


elo  z  *  2(Fa-3)fa 


+  terms  of  order  >  t 


If  we  introduce  coefficients  of  variation  by 


these  fall  into  the  form  quoted  in  Section  11. 
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LET AILS  OF  SECOND  EXAMPLE 


27 .  The  responae 

We  next  consider  a  symmetrical  ir-section  attenuating 
network  operating  between  equal  image  impedances  R.  The  attenua¬ 
tion  can  be  shown  to  be  given  by 

z  - 1  +  m  +  2  *b  !r^  +  ^  j+ 1  ^  + 


where  Rg  is  the  resistance  of  the  series  arm,  and  R^  and  Rc  are 
those  of  the  two  shunt  arms.  The  ideal  values  of  Rg  and  RQ  are 


RA  -  Rc  -  R  , 


a2-l 


where  a  is  the  design  attenuation, 


If  we  place 


ra  =  O  +v 


”b  ’  R~4r  O  +  V 


Rc  -  R  -I!-v  t1  +  W„) 


the  attenuation  becomes 


z  =  a  {1+wArl  (1+WC)’1  1  +  sir  +  <V*C> 


+  WB  <  +  WAWC  ~ Jr1  +  WB 


and  we  are  almost  ready  to  find  derivatives. 


a*  m  mmm .=■ 

' ■*  •  ..  ,  ’  ,  . 


ii£, 


SjHpji 


W»*^!*S«isefegabi*6*^ 


-  vil-k 


* 


Accurate  di f ferentiatlon  will  however  be  swifter  and 
easier  If  we  rewrite  this  expression  as  a  sum  of  products  of 
functions  of  the  single  W?s,  namely  as: 


3  1+WA  l4-WC  +  2  a+JL  1+VA 


— z t  /  1  ,  _  a+3 _ A  1 

i+WA  ’  B  1 +WQ  a  2(aJl)l+WA  T+VT“ 


I  _ _C  u-1 
1+WA  1+W„  +  2 


1+WA  "b  1+WC 


2  1+W, 


w  -T.-C  +  af+4a-l  WA 

B  1+W 


+  Vi  _!a_  w  wc 

T+W^  WB  T+W  ~ 


and  note  that  the  values  of  the  successive  derivatives,  at  W  =  0, 
of  (l+W)  and  W/(l+W)  are  as  follows: 


order  of  diff»n: 


(l+W)-1 

W/(l+W) 


We  thus  obtain,  writing 


z  =  "<V  V  V  . 


and  introducing 


\  “  ' , 
!i 


'Sfewftsi'S  ■  .  . 


Vii-3 


a  3  3-2 

P  =  1?  a+T~  •  y  ■  -4- 


cuo+lj  ' 


which  satisfy  20?  a  62 


HA  =  H  >  -g,  H..  -  Hrr  =  2(3,  H 


AA  *  »CC  =  HAAA  =  HCCC 


haaaa  -  Kcccc  =  24f)-  haaaaa  =  Hccccc  =  -12°P 
HAC  =  7’  KAAC  "  hacc  =  •2’,<  haa an  ’  V,  H..., 


:aacc  *  haaac  =  haccc 


haaacc  “  '12/< 


HAB  “bc  -  -3’  haab  “  HgCC  =  26,  H 


AAAB  ^CCC  “  "66> 


^BB  ABB  ^AABB  ^EBC  ^BBCC  ~  ^ 


ABBC 


We  are  now  almost  ready  to  write  down 


formulas 


the  generalized  propagation 


It  is  convenient  to  first  evaluate 


with  the  results  given  in  Table  8 


certain  coefficients. 


L*j|  t— 1 
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TABLE  3 

Values  of  Coefficients  In  the  Second  Example 

Value  When 


Coefficient 

ab=AB  or 

CB  ab=BA  or  BC  ab=ACorCA 

HaHabb+Kab  +HaabHb 

2HaHaabb+5HaaHabb 

62+4(35 

62+436 

y2+U(3y 

^"aab+Haab1^ 

-262-4B6 

0 

-6fiy-2y2 

6Wab 

12B26 

12P26 

6fi2y 

3HaHab2 

2 

-336^ 

6B62 

2 

-337 

lHa2«abb+3HaaHabHb+3HaHaabHb 

-24p26 

12£26 

-15 P2y 

Value  When 

A 

a=B 

a_ 

=A  or  C 

Ha2 

43  2 

p2 

Vaa  "  T5HaHaaaa 

0 

-2B2 

-iff  H  =  _  iu  H 

3  a  aaa  3  aa  aaa 

0 

232 

_£h  d 

4  aa 

0 

q2 

3 

Ha3 

83  3 

-  33 

KXa 

C 

333 

2Ha2,1aaa  '  lHaHaa2 

0 

-333 

Jslng  these  results }  the  ■fo’pnrulas  become 


ave  z 


var  z  = 


«  a 


+  3  (a/  +  ac2) 
-  p  (^“a3  + 


B  C 

4  .  ..  4, 


+  P  (rA^A  +  rcac4) 


2  2 

+  ?aA  aC 


'  P  (GActA5  +  an°r5) 


cwc 


2_  2 


+  7  (rA oA  +  7cac)a/oc 


+  .  .  . 


e2  (c/  +  4cb2  +  oc2) 

-  2P2  (?A°A3  +  ^0ac3) 

+  2p2  <rA°A4  +  Vc4) 


+  P‘ 


(rA-1)oA  4  (V^C4) 


+  (62  +  4P6)ob2  (aA2  +  Op2)  +  (72  +  4g-,)oA2 

-  2p2  <Va*  +  °C"C5) 


2£‘ 


(aA~yA^aA5  +  (ac~^aC5 


(SB2  +  4P6 )(rAaA3  +  rcac3)oB2-(2y2-6py)(7Ao 


+ 


vxi-6 


ske  z  *  ~  8yBaB3  + 


+  33: 


+  I2p5  ( <ja  +  oc2)gb2  +  6P2yGA2°c2 

r  -| 

-  3e3  (oA-rA)oA5  +  (oc-rc)oc5  i 

.L  — J 


-  33' 


{Gh’2yk)ah  +  (Gc'2^c)ac 


-  33  6  ‘ 


(’'A- 1 )  ''a3  CB2-2  '^B- 1  >  °B3 1  °A2+ae2 }  +  ( >  C-1  JaC3crB 


7| 


2  2 

0C 


-  3@r2  (jv1*  V(V1)oc]  <V 

-  6P2&  f(^A'’A3+i*709c3)aB2+27BaB3(OA2+Oc2)j 


J 


2_  2 


-  13P  °c 


I  •  a  • 


elo  z  =  0 


^rA"3^0A4  +  i6  (rP“3)  Cg  +  (rn-3) 


-  43 


(aA-4>A)oA5  +  ^C’^c'V 


+  .  .  . 


To  bring  these  results  into  a  more  perspicuous  form, 
we  change  to  a  more  natural  scale  by  putting 


’*’**ww» WMWHMWlcaW  Wte 


B8S  Hjggj 


'<}$,  . 

mu  $i  ■ 


•v.: 


K  ■- 


TO&\1y 


'  „;S.- 
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WA  *  0,1  V  WB  *  C‘a  V  WC  ”  °‘1  VC 


z  =  H(Wa,  Wb,  Wc)  *  h(wA,  WB,  wc) 


so  that  the  standard  deviations  of  the  w»s  will  be  of  the  order 
of  unity  or  somewhat  smaller.  Notice  that  0,  y  and  5  are  sub¬ 
stantial  fractions  of  a,  for  example 


a: 

1 

2 

5 

10 

uo 

3: 

0 

0.17a 

0.33a 

O.^la 

0.50a 

y: 

C 

0.02a 

0.11a 

0.17a 

0.25a 

6: 

0 

0 . 08a 

0 . 27a 

0.37a 

0.50a 

so  that  we  may,  roughly,  think  of  a,  0,  y  and  6  as  of  about  the 
same  order.  The  resulting  formulas  were  already  given  in  Section  12 
28 •  Transforming  Individuals , 

If  we  revert  to  the  original  expression  for  the  trans¬ 
mission 


1  +  2R  "b  +  ^  +  +  I  +  J 


'_L  ,1  .  \ 

*A  RC  W] 


and  consider  the  effect  of  varying  RA>  Rg  or  Rc  alone,  we  see  that 
z  is  linear  in  1/RC.  Thus  the  appropriate  transforms  for  R  and 

f* 

Rc  arise  by  replacing  resistances  by  conductances,  while  the  ap¬ 
propriate  transform  for  Rg  is  the  identity. 


If  we  place 


Ra  =  P  “a?T"  '2  +  V 


%  R  ”2q2  ^  +  vs) 


1  a-1 


Rc  R  a-1-! 


(J-  +  V_) 


then  the  attenuation  becomes 


2  =  a 


+  2  Sf  <  V2VV  +  T^THriT  Vc  +  jfsir 


1 2 


~f  VB 


4- 


y  ^ 

VAVBVC 


orJ  ln  abbreviated  notation 


a  +  P(Va+2Vb+Vc)  +  rVAVc+6VB(VA+Vc)  +  yvAvBvc 


and  if  we  write  z  =  f(vA,  v0,  V )  then 


we  have 


fA  *  fC  “  fAC  ~  y 


=  23,  f 


AB  BC 


fABC  ~  7 


while  all  other  derivatives  vanish  for  V.  =  VQ  =  y  =  o.  Thus  the 

ABC 

generalized  propagation  formulas  reduce  to  : 


ave  z  =  a 


var  z  -  p2(oA2  +  4ob2  +  0(?) 

,  .2,  2  .  2 


3ke  z 


+  -  °CW  +  ?ws 

+  .  ,  , 

"  *3(  W  +  8W  + 

+  12326(aA2  +  crc2)aB2  +  6327c?A2ac 
+  3362  (7A“l)cA8  +  (7c-l)ac8)  0] 

Lj~  — j  _ 

+  33  72  +  (7c-l)oc  cA2a( 

+  6^2(7B“i)aB3(aA2  +  ac2)  ~ 


2„  2 

A  aC 


4* 

*  *  *  # 


lo  z 


e4  (rA-3)aA4  +  i6(rB-3)<JB4  +  (rc-3)oc 


With  the  possible  expectlon  of  the  formula  for  ske  z,  where  the 
reduction  in  number  of  terms  is,  though  still  large  in  number  only 
about  half  of  the  number  present,  the  gain  over  each  of  the  previou 
formulas  is  notable. 

Rescaling  again  by  writing 

VA  =  °-1  VA 

VB  =  0.1  vB 

Vc=°.lvc 

we  obtain  the  formulas  already  given  in  Section  12. 


*£*•=*- 
¥'  h&  J».  *§ 


DFTAI LS  OF  TRANSFORMATION  OP  RESPONSE 


As  Indicated  in  Part  II,  the  attack  by  transformation  can 
involve  both  transformation  of  the  individual  [component!  variables, 
which  is  naturally  directed  to  annulling  the  higher  unmixed  deriv¬ 
ative  values,  and  transformation  of  the  response  [system]  variable, 
which  might  naturally  be  directed  at  either  increasing  the  normality 
of  the  response  distribution  or  at  reducing  the  size  of  some  of  the 
terms  involving  mixed  derivative  values.  It  appears  that  trans¬ 
formation  of  the  individual  variables  is  not  only  easier  but  rathor 
more  efficacious  Thus  we  will  be  wise  to  plan  to  transform  the 
individual  variables  first.  This  will  be  true,  even  although  an 
incautious  transformation  of  the  response  would  disturb  the  desir¬ 
able  situation  obtained  by  transforming  individuals,  since  we  can 
be  cautious,  and  arrange  for  a  compensating  transformation  of  the 
individual  variables  to  accompany  the  transformation  of  the  response. 

We  thus  contemplate  the  following  sequence  cf  situations; 


(z  h ( w ^ *  Wgj  w^.) 

(perfectly  general 
^Z  =  g(v^;  Vgj  . . * »  v^) 

(v  a  function  of  w 
'  a  a 

iSaa  =  gaaa  ”  . . .  0 


-■  •-  .  -  --  . 

■*  >%  .  fa#  ’  ' 
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jy  *  f(ur  u, 
<y  SI  ?(«) 


va  '  va(V 


k 


aa 


'aaa 


0 


Here  we  may  as  well  require  that  the  first  derivative  values  of 
ee  h  v_  with  respect  to  w  .  of  y  with  respect  to  z,  and  of  each 
v  with  respect  to  u^  all  be  unity  in  order  to  keep  all  leading 
terms  the  same.  We  naturally  denote  the  latter  two  derivative 
values,  evaluated  at  the  corresponding  average  points,  by  <p 1 
and  v  1 ,  We  may  do  this  because  we  are  dealing  with  functions 

d 

of  single  variables.  We  shall  use  cp",  cp  * 1  ' ,  cpiv,  V  ",  V  *  »»,  v  1 

B.  a  ci 

for  the  corresponding  higher  derivatives. 

We  have  already  discussed  the  first  step  in  some  detail 
(Section  6ff.)  so  that  we  need  here  to  be  concerned  with  the  sec¬ 
ond  step.  If  we  are  to  carry  out  the  second  step  in  practice,  we 
need  to  be  able  to  find  out  about  vQ(u  )  and  about  f,  *s  in  terms 
of  cp  (presumably  in  terms  of  the  derivative  values  of  cp)  and  the 
sab*  Tills  we  shall  do  next. 

30.  Transfer  formula  details 

Suppose  now  that,  starting  with 


z  =  g(v1#  v2,  . . vk) 


s  m  s@s**a*>»!* 


mm 


-  Vlll-3  - 

where  saa  =  ^aaa  =  saaaa  =  =  we  t rans f orm  the  response  by 

y  ~  <p(z ) ,  and  wish  to  ms.ke  transformations  on  the  individual  vari¬ 
ables  va  =  va(ua)  so  that  in 


y  -  <p(g(v1(u1),  v2(u2),  ...  vk(uk)))  =  r(u1}  u2 ,  ...  u  ) 


vj  ■  will  have  f&a  ~  faaaa  -  -  0.  We  surely  need  to  have 

formulas  expressing  the  (successive  and)  cross  derivative  values  of 

f  in  terms  of  the  successive  and  cross  derivative  values  of  g  and 

the  derivative  values  of  <p.  Information  on  the  functions  v  (u  ) 

a'  a 

will  also  be  helpful. 

There  is  no  loss  of  generality  in  assuming  that  the  first 
derivative  values  <p',  v^,  v2‘,  vk  *  are  all  unity.  With  this 

convention,  the  higher  simple  successive  derivatives  of  f  are* 

fa  “  *'  Va'  “  Sa 


f»a  '  9"(gV  ')"  +  <P*  g„„(v  ')d  +  o' 


aa '  a 


9 1  Sava  "  =  9 "  ga  +  Saa  +  gava  11 


f  =  q>  1  (g  v 
aaa  y  ve>a  a 


'"(gaV)3  +  39"  eaaga(va')3  +  39"  gava'va" 


+  <*'  V(v»')3  +  39’  Saava"va'  +  9'  g0vo" 


9'"  g3  +  39”  SaaSa  +  39"  gfv„"  +  g. 


+  3  ®aava "  +  Va'" 
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SHciLC. 


$  (%  V  * )  ’  +■  6m  ’*  a  etC t \*  »  V*  a  #5- »  »  „-J„  "7 ,,  I  \ 

u  n  ’  v’aa,Ja  ‘  *  a  4  ",'r  &a¥a  / 


a  a  a 


'  saaasi. ( va  ‘  >  +  18<P  "  <5„  „S„  ( v.  1 ) 


aa“a '  a 


+  3q5"  safvaH)~  "  4(P  ’  «av-  ” 


"■  V  ’  V  "  !  —  m  1  tr  (  v  1  ) 

a  a  a  4  --aa a  a  a 


+  V  ®aaa^va ' ^~va "  +  3* '  «aa(va"^  +  4*'  S^VV 


'aa  a  a 


+  <P'SaV  +  3*"  Saa<V> 


1>1V  Sa  +  V"  Saa«a  +  6tP'"  Sava"  +  V  Sat.aSa 


+  "  ®aasava "  +  3<p"  Sa(va")2  +  V  &J" 


:aaaa  +  5  saaava  +  3  &aa^va"^  +  4  saava"‘ 


iv  2 

+  g  v  +  3g>  g 
a  a,  °aa 


If  now  we  set  these  successively  equal  to  zero,  starting  with 
second  and  using  g?a  =  gaaa  -  sa£aa  -  =  0,  we  obtain 


0  =  f"  «a  +  ®aVa " 


0  =  g"'  g3  +  3<p "  gfv  "  +  g  v  " * 


0  -  q>lv  g4  +  6<p"'  g3va"  +  3<p"g^(va")2  +  4ip"  g2va"-  +  gavalv 


and  similarly 


yr^T-r, 


~  ©*  S'  +  10®*^  (E^v  "  4-  lOfD  *'  *  M  *  4,  C,  n  T~  Iv  ,  _  ,,  V 

v  °a  T  ®a  a  ^  ®*,*a  +  3<P  2pv«  +  S  v 

a  a,  cl  d 


+  15<pf,»  g?(v  ")2  +  10®"  g2v  ** * v  » 
d  a  a  a  a 


whence 


V  “  -(sa<P")> 


"  I 


S2a  -  3?"  Sava" 


■  -<e%"’)  +  3(g.tp")2 


=  gr(V"  +  39V), 


v. 


lv  -3_iv  *  2  ill..  „ 


a 


=  -g^±v  -  6  g^“va"  -  3<p"  sa(va")2  -  49"  g„v= 


'a  a 


*  -<^1V)  +  10(g2,,"')(ga9")  -  15(ga9")3 


=  +  10qj"'9"  -  15(p VV). 


and 


v  4  r  V  1  o  1 

va  =  ga  -CP  +  l5cP  (P"  +  10cp  11  fqp " 1  -  105qp "  1  ( <p") 2  +  105(q>")^]. 


with  these  formulas  In  hand,  we  are  prepared  for  the  cross  deriv¬ 
atives,  and  find 


I 

i 
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"  V1  8avaVb'  +  S»kv»'v- 


ab  a  b 


^  sagb  ' 


fabb  “  "'"  Sava'(sbvb')2  +  2<f"  gbgabV(V>'  +  <p"  Sava  'gbb  <vb 


+  9"sava'gbvb"  +  f'  5fl^va  1  f  VK  * ) 


abb  a  1 'b 


+  qp 1  g  .  v  1  v,  " 
y  6ab  a  b 


■  gagb  +  2’>"  ®abgb  +  ?"  ®agbb  +  ®agbvb"  +  Sabb^abV 
"  ®abb  +  *”  gab®b  +  ?"  ®agbb  +  gagb^'"  ‘  <?>") 


fabbb  "  '*'1V(eava,)(gbvb')3  +  V  gava'gbgbb(vb'i3 


+  39"'  SaVgbvb'vb"  + 


+  <?"  S„Vghhh(vh')3  +  39"  Sava'8bb(vb,)(vb") 


+  39"'  Kabgbva'(vb')3  +  39"  S^g^VK')3 


abb£b  a  v vb 


+  V  gabgbbva'(vb')3  +  6?"  S=hS*WV 


absb  a  ”b  vb 


+  9"  Sava'gbvb'" 


(formula  continues) 


mi -7  - 


+  ©  * 


(vQ*)'vvn  3 


=*abbbv  a  ;  b 


+  39 1  g  v,wv  'v,  *v.  *'  + 
°abb  abb 


+  *'  gabva  *vb 


91V  gag3  +  3<p" 1  gaSbSbb  +  3cp!''  gQg^vb"  +  icp"1 


’a°b  b 


5ab°b 


+  *"  gagbV  +  SaSbbb  +  3? M  gagbbvb"  +  39"  gabbgb 


+  39"  gQ*g 


ab&bb 


+  6*"  *ab*bvb"  +  gabbb  +  3gabbvb"  +  gabV" 


gabbb  +  3<P  gabgbb  +  ^  +  (29  '  -  39  9  )gav,g 


a°bbb 


Jab&b 


+  3(9"'  *  9  "9  ")gaSb8bb  +  (9  V  "  ^9  "9 " 1  +  39 '9  !>9")gagb 


’abbb 


+  (2?'"  -  3<pV)gabgb  +  (qrv  -  VV  +  3?VV')gag2 


mmmtmm gftghtfet 
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1 


4  gabgbv/  -  gaagbb  +  2f”  gs«2fcb  +  <?’* 

+  *"  gaagbV  +  2<PH  gasabV  +  *"  gagbW 

H  gabbva  +  gaabvb  +  gabva  vb 


tr  jT  v  n 

®agbbva 

+  gaabb 


gaabb  +  ^  gaabgb  +  ^  gagabb  +  ^  gaagbb 

+  +  (9"’  -  *V)gfgbb  +  (« P"'  '  <p'V')saag^ 

+  (V1'  -  3<p"<p")ga£absb  +  (<plv  -  2<p",<p"+(P’VV')gagb 

‘  gaabb  +  9"<8aal.8b  +  gagabb>  +  V  «ab 

+  <V"  -  3<pV)  gagabgb  +  (tp1V  -  2?"’<p"+9,VV’)g2g2 

(In  each  case,  the  first  equals  sign  is  followed  by  the  general 
expression,  the  second  by  the  result  of  using  unity  for  tp 1 ,  v  1 
anc*  vb  *  *  third  by  the  result  of  further  substituting  for 
va">  vb">  and  vb"S  and  the  last  by  the  result  of  using 

gbb  =  gbbb  =  0  * ) 

These  results  are  to  be  supplemented,  when  needed,  by 

1  P 

va  «  ave  va  +  ua  +  5  Sa(V)u‘ 

+  5  ga(-<p"'  +  3<p'V) 

+  |4  ga(-'plv  +  10<p"'  <p"  -  15<pVV') 

+  .  .  . 


i  ■  V  ■-& 

,y“ '■ ‘v;'-  ■ ,j '  -*  ■' 
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which  Inverts  Into  the  Taylor  expansion 
ua  -  (va  -  ave  va) 

+  i  ga  " f va  "  ave  va^? 

1  2  it  i  /  \  3 

+  B  ea  *  (va  “  ave  va' 

+  k*l  -  ave  va)4 

+  .  .  . 

as  it  obviously  should,  since 

ua  Mg(ave  vy  ave  Vg,  .  ..,  va.  . ..,  ave  vkJ) 

-cp(g(ave  vv  ave  Vg,  . ave  v&,  ...  ave  vk))]. 

When  the  formulas  giving  f's  in  terms  of  g's  are  trans 
lated  into  formulas  giving  s's  in  terms  of  t's,  we  obtain  the 
formulas  already  given  in  Section  14. 

31.  Improving  normality 

We  are  now  ready  to  consider  the  possible  modes  of  use 
of  transformations  of  response.  First,  consider  the  possibility 
of  choosing  q>  so  that  the  normality  of  g  *  <p(z)  will  be  improved. 


<*# 
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The 


leading  terms  in  ske  y,  where  y  =  f  (u-, ,  u.  ,  . uK)  are 


Z  yr\  +  Z  (6s  ,  )t  t 
' a  a  ab'  a  b 


Jr 


=  X  7ar|  +  X*(6tab  +  6v")r'l4 

where  the  t  are  calculated  for  the  u  .  (Hopefully,  but  not 

cl  3. 

certainly,  the  next  term  will  not  be  Important.) 

In  order  to  reduce  the  size  of  this  approximation  to  the 
skewness,  we  may  proceed  along  two  different  paths.  We  may  try  to 
make  the  effects  of  some  terms  compensate  others,  or  we  may  try  to 
make  as  many  individual  terms  small  as  we  can. 

The  difficulty  in  seeking  reduction  by  compensation 
stems  from  our  lack  of  knowledge  of  the  a  and  y  .  Setting  the 

cl  3 

expression  equal  to  zero,  and  solving  for  <p"  yields 


?  r  *  2  ? 

<  +  6z  t  VT“< 

_n  _  'a  a  abab 

cp  =  .  — — r-g— 3 - 

6z  Ta  Tb 


which  is  rather  unmanageable  if  the  Ta  and  y are  either  not  pre¬ 
cisely  known  or  may  vary.  In  some  circumstances  such  reduction 
might  be  practical,  but  its  general  utility  is  no  better  than 
doubtful , 
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The  difficulty  in  seeking  reduction  by  reducing  in¬ 
dividual  terms  Is  two-fold.  There  are  terms  which  remain  unaffected. 
And  not  all  affected  terms  vanish  for  the  same  value  of  <p".  The 
straightforward  condition  for  eliminating  a  particular  term  is, 
however. 


<P 


It  ~ 

<V 


Clearly,  most  attention  should  be  given  to  large  values  of  TaTfe. 
If  we  weight  by  we  find 


2  Ta 

- ? 

Ta  Tb 


2 

Tb 


We  may  not  know  the  t*s  well  enough  to  use  them  at  all  accurately, 
but  we  may  get  a  rough  idea  of  a  possibly  useful  cp"  from  this 
relation. 

Further  arguments  against  the  usefulness  of  attempts  to 

improve  normality  can  be  based  on  the  relative  size  of  the  y&  and 

the  x  .  If  the  former  are  much  the  larger,  then  cp"  will  be  large, 
a 

and  other  terms  in  ske  y  may  well  be  important.  If  the  are 
large,  the  same  may  be  said.  If  all  are  small,  and  the  y &  are 
smaller,  then  the  skewness  is  small  —  why  do  we  concern  ourselves 
with  it? 


**»*•«! 
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32,  sm 


Clearly  we  must  choose  which  formula  is  to  be  simplified 
in  which  way,  (Obviously  we  are  going  to  try  to  simplify  the 
variance  formula  —  but  it  may  be  worthwhile  to  rationalize  this 
choice.)  The  arguments  of  the  last  section  speak  strongly  against 
simplifying  the  formulas  for  skewness  or  elongation.  The  second 
term  in  the  formula  for  the  average  involves  values  of  4^  deriva¬ 
tives  f  faabb J  •  So  our  attention  is  directed  to  the  variance, 

where  the  terms  next  after  the  first  involve  t2t?  .  If  we  look  for 

a  b 

only  the  lowest  order  of  differentiation  we  are  led  to 


^  2  2  2 
2  sab  TaTb 


and  matters  proceed  as  in  the  second  approach  to  increased  normality, 
since  larger  values  of  t2t^  deserve  greater  weight.  If,  again,  we 
can  disregard  variations  in  the  sizes  of  the  a  ,  weighting  by 

EL 

2  2 

sagbJ  which  leads  again  to 


Z  gagbgab 
— 

S  ga  sb 


may  be  of  use.  All  this  concerns  only  part  of  the  coefficient  of 
2  2 

TaTb'  however* 


% 
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If  we  concern  ourselves  with  the  whole  coefficient  of 

2  2 

*r&  Tb  '*  we  face 


2 

S  ,  +  S  v.  "t*  S 

aab  ab  abb 


which  is 

Ik 

'l 


I 


4 


^aab  +  ^  +  ‘abb’  +  4tab^"  +  (V"  '  <p"<p") 

By  choosing  <p 11  *  and  9"  appropriately  we  can  make  this  vanish  for  any 
two  pairs  (a,b)  or  for  any  two  sets  of  equivalent  pairs.  As  the 
examples  will  show,  this  may  sometimes  be  useful. 

33.  Power  transformations 

The  natural  transformations  to  consider  are  often  of  the 

form 

y  =  A(z  +  c)p 

whence 


<P*  =  pA(zo+c)p_1  ^ 

cp"  =  p(p-l)  A(zo+c)P"2, 

l. 

cp"*  =  pfp-l)(p-2)A(zo+c)p'3, 


q>  ’  *  (p-i)/(zr.+c) 


(p-l)(p-2)/(z  +c)‘ 


2qp"»  -  cp "9 "  =  [2(p-l)(p-2)  -  (p-1)  ]/(zo+c) 


=  (p-l)(p-3)/(z  +c)‘ 


whence  p  may  be  found  from 


1  -  P  = 


11  l  MS 

cp  1  -  cp  qp 


z„  +  c  as 
o 


+  c  =  - 

O  -cp 


and  A  as 


(z0+°) 


WSWiMwpffto  mu; 


VTII-16 


to 

* 


It  xs  worthy  of  remark  that  the  case  p  *  0  corresponds 


y  -  (z0+c)  log(z+c) 


with  <p>  =  1,  qp"  =  -l/(zQ+c),  <p"»  =  2/(zo+c}2,  and  <piv  =  -6/(zo+c}3, 
as  might  have  been  expected. 

Similarly,  as  p  *4±  M)  the  transform  approaches 


(  i.  ) 

z0+c 


Z  +C 
o 


■)(z-20) 


* 

where 


* 


Especially  since  the  values  of  p,  c  and  A  may  be  quite 
sensitive  to  the  values  of  cp"  and  cp"»,  it  is  useful  to  have  formu¬ 
las  for  qplv  and  cpv  in  terms  of  cp"  and  cp1.  Both  for  the  general 
and  limiting  cases  of  power  transformations,  we  have 


* 


so  that  the  more  complicated  formulas  connecting  s's  and  t’s 
become 

sabbb  =  fcabbb  +  (2<p”  '“3<p"(p")tab  +  (2q>"  f~qp ”cd  -  <p  "(4<p "  »-3<p "<p  ") 

saabb  =  ‘aabb  +  «’"(taab+2tab+tabb)  + 

+  (2<p"'.»'V')(-^-  -  ,») 

saaabb  =  fcaaabb  +  6<p"  tab  taab  +  V  ' ' -<P "<P ") taato 

+  (29’'<-3<pV)tabb  +  6(9".-(p'V’)ta2 

+  |l2(<p'")2  -  22q>"VV’  +  9(<p'V')2]  tab 

+  l6^”’)3  *  16(<P"’)2  9'V  +  13*"*(q>V)2-3(q>V)3j 


mmmm 


We  are  really  more  concerned  with  formulas  for  the 
coefficients  in  the  propagation  formulas.  Combination  of  the 
formulas  just  given,  and  the  use  of  suitable  abbreviations,  lead 
to  the  formulas  already  given  in  Table  5  (Section  18). 

If  we  seek  to  fix  <p"  t  nd  qp "  *  by  arranging  for  the 
individual  coefficients  of  and  to  both  vanish,  then, 

solving  the  simultaneous  linear  equations,  we  must  take 


where 


qab  taab  +  tab  +  tabb  ' 
qcd  ~  tccd  +  ^cd  +  fccdd  ’ 


♦ 
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TRANS  FT  HAT] ON  OF  RESPONSE  IN  THE  £XAKPLES 

3^ *  Petal Ic  for  first  example 

After  we  had  made  the  Individual  transformations, 

we  had 


S 


21-1  s2i 


g 


21-1,21  *  ZRJ 


with  all  other  higher  derivative  values  zero.  Hence 
t2i-l  21  =  while  all  other  t's  vanish. 

If  we  try  to  reduce  Individual  coefficients  in  the 
skewness,  we  have,  when  we  equate  <p"  to  -t  ,  for  various 

Q.  D 

choices  of  (a,b) 

<p"  %  -0.1  (j  occurrences) 


<p"  ~  0.0  (2j(j-l)  occurrences) 


and  clearly  we  are  likely  to  orefer  to  keep  cp"  =  0  unless 
j  =  1.  The  weighted  solution  gives 


_n  'v»  1 

*  "  “  10(2 j -1 ) 


which  is  clearly  quite.1  near  zero. 


; 
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9 

For  the  most  plausible  approach  we  find 


* 

0  =  0.01  +  0, 4<p'  +  2q?H  -  <p'V'  { j  times) 


0  -  0  +  0  +  2qp"  1  -  q>V'  (2J(J-1)  times) 

and  we  can  satisfy  both,  if 


<?"  =  -1/40 
<p"'  =  1/3200 


* 


which  corresponds  to  a  power  transformation  with  p  =  3, 
2Q  +  c  «  -80,  A  *  l/(3(-80)2),  namely  to 


y 


(z-10,1-80)3 

19,200 


for  which  <p'  .  l,<p"— 1/40,  cp',f,=  1/3200,  and  q>lv  .  <pv 


»  0. 


The  resulting  values  of  s^,  etc.  are,  when  a  and  b 
are  paired  (that  is  when  (a,b)  *  (2i-l,  2i)  or  (2i,2i-l)), 


•* 


*► 


£:  mmm)$ 


t 


s  ,  » 

ab 

^aab 

3abbb  = 

saabb 

s  ,  .  - 

aaaDD 


TV.t 


c.l  -  0. 025  *  0. 0i  5 

3abb  =  -0.025(0.1)  +  =  -0. 00262 + 

^aaab  =  ^3200  "  I^OO^0,1^  +  (4 0)3200 

+  __ - 2__  =  -  0.00014" 

(-40)3 

-0. 025(2) (0.1 ) ^  +  (3±_  -3^) (0.01)  -  TTCJ§5ro 


(-40) 3 

-0.00056+ 


saabbb  6 ^3200  “  IHJo^0*1)2  +  (’fifo)  (3200) 


I 


4  i_)(0.l) - i *■  + - - - ^_T_ 

(40)-3  (3200)*  (3200)  (40)2  (4or 


A. 


4  0.00009 


These  lead  to  the  following  coefficients  (still  to  5  decimals) 


k  saabb  “  '  jg  saaabb  -  +  0. 00001 + 

e 

2 

3aab  +  3ab  +  3abb  "  0 

?  3aabb  +  sabsaab  +  J  saaab  "  -°-000^ 

6sab  ’  °'45  3sab2  -  .01688' 

3saab  +  2  3abb  '  -°-01265' 

and  hence  to  the  formulas 
ave  y  *  -80 


-0.00014  z  t21.iT£12 


2  2 


+0.00001  -t-.  +T-  )r  *•  — 

w2i-lT2i-l+72in21~l  T2i 


var  y  «  St 


+  terms  of  order  >t  , 
2 


2_  2 


•0.00052  Zb'gi.jTg^  +  ^21T2i  ^t21-1^t21 


+  terms  of  order  >tG  , 


ix -5  - 


ske  y 


elo  y 


ponding 


ave  z  = 


var  z  = 


S?aTa 


+  0.45  2  ^21-1  T2i 


+  0.01688  2[(721_1-1)t21„1+(721-1)t21]t21_1  t2i2J 


2  2 


0.01265  ^721-1^21-1  +  y2iT2i^T2i-l  T2i 


+  terns  of  order  >  t  , 


+  terns  of  order  > 

These  formulas  are  to  be  compared  with  the  corres- 
formulas  before  transformation,  which  are 


lOj  +  no  other  terns. 


+  °-01  ST21-I%i2 

+  no  other  terms. 


'  %* 
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SJ€2  Z  *  Zy  x  *' 

'a  a 


+0.6  Zto,  A  2 


k2I-I  ‘21 


+  °-03  2^2i-AAi-i  -  Ai^Ai-iS*2 


+  terms  of  order  >  t  , 


elo  z  -  2(ra-3)ra 


+  terms  of  order  >  A 


These  two  sets  of  formulas  are  to  be  assessed  in  view  of  the 
fact  that  ira  =  ~  oa,  so  that  t's  as  large  as  1.0  are  possible 
but  unlikely.  Values  between  0.5  and  0.05  are  more  plausible, 


Thus 


0.0005  +  y2i,21U2^\  2 


is  not  likely  to  be  more  than  a  twentieth  the  magnitude 


2  2 

+0.01t21_1  t21  and  might  be  considerably  less. 


In  general  terms,  the  transf ormation  has 
(i)  added  some  minor,  almost  surely  negligible,  terms 
to  the  formula  for  the  average. 


ifjrlL . 

4?f  :  -  -%S'  1 


.  ^ 


ftsMwK*  m  :;T 


(il)  reduced  the  most  important  correction  term  in  the 

variance  from  perhaps  1%  cf  the  leading  term  tc  per¬ 
haps  0.05^  -  from  almost  certain  negligibility  to 

utter  unimportance. 

(ill)  modified  the  correction  terms  in  the  skewness,  and 

has  had  some  slight  tendency  to  reduction  of  coef¬ 
ficients. 

Overall,  transformation  has  done  what  it  was  asked  to  do  —  it 

has  eliminated  the  most  Important  correction  term  in  the  formula 

for  the  variance.  However,  in  this  instance,  the  correction 

term  in  the  variance  was  not  important  enough  to  warrant  much 
effort 

35 •  A  Salutary  Example 

In  searching  for  an  example  showing  how  much  trans¬ 
formation  of  the  response  could  matter  if  „e  had  picked  the 
wrong  terms  of  response  initially,  we  may  well  ask  what  would 
have  happened  it  we  had  started  with  a  response  proportional 
to  the  cube  of  the  delay  time  of  the  filter.  We  can  easily 
calculate  the  reduced  derivative  values  by  taking 


300j' 


when  =  1,  „»  =  !/5j,  ,  1/50j2>  9iv  _  =  _  _  _ 


The  sab  et0  >  which  now  do  not  refer  to  well  chosen 


terms  of  response,  become 


sab  »  0.1  +  0.2/j 

aat  «*  0.02/j  +  0.02/j2  -  0.04/j 2  •  0.02/j  -  0.02/j2 

abbb  =  (0.04/j2  -  Q.12/j2)Q.i  -  4(0.004/j3  +  3(0.008/j3) 

»  -  0.008/J2  +  0.008/J3 

saabb  “  (°-2/J ) (0.02 )  h  (4(G.02/j2)  -  3(0.04/j2))  (o.l) 

-  2(0.004/j3)  +  (0.008/j3) 

=  0.004/j  -  0.0 04 /j 2 

The  leading  and  first  correction  terns  in  the  variance 

become 

var  y  =  Zt  2 

ci 

+  (0.01  +  0.03/j)  2  ^ 

4~ 

«  «  •  » 

For  very  small  J  (namely  J  =  1 ),  beginning  with  the  cube  of  the 
delay  time  as  the  response  might  make  the  first  correction  term 
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relevant,  but  for  larger  j  (namely  j  >  2}  the  first  correction 
la  small  enough  to  be  unimportant ,  even  when  this  very  poor 
choice  of  terms  is  made. 

36.  Details  for  second  example 

After  we  had  made  the  individual  transformations, 

we  had 


with  all  other  t's  vanishing. 

If  we  try  oo  reduce  individual  coefficients  In  the 
skewness,  we  have 


If  'V 

■  ^5  -v 


1 

a 


a 


„  ti  ~ 

-cp  ~ 


l+l 

a 

a 


( one  e ) 

(twice ) 


and  the  choice  a)  «  -i/a  seems  likely  to  be  a  helpful  compromise. 
The  choice  p  *  0,  c  =  0,  since  zQ  =  a,  will  do  this.  Thus 


y  *  a  log  z 


a  .  ,  ?■ 


%  r  : 


,  %  % 


,,7~- 
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* 


* 


f 


* 


t 


is  to  be  considers  .  The  weighted  solution  differs  from 
<p”  *  -1/a  by  a  small,  presumably  negligible  result.  If  we 
try  to  delete  mixed  correction  terms  In  the  variance,  we  find 


I  +  1 

a  + 
l  a 


+  4 


Si  i)  „ 


-+- 

a  ) 


If  _  If 

-cp  cp 


0  (twice) 


>1 

>a 


+  4  (± 
)a 


— 5-icp"  +  2cp,TT  -<p"<p"  *»  0  (once) 
a-5  ) 


which  leads  to 


cp  =  - 


+ 


2a  4aS  '  4? 


-It  '  _  0 

<p  -  -T5 


11  3  +  7  +  1 

8a‘  32a  32a4  16c?  32a6 


We  might  reasonably  hope  to  obtain  good  results  with 
cp'  =  -l/2a  and  qp"  =  5/8a  which  correspond  to  a  power  trans¬ 
formation  with 
A  =  3(W3)2/3.  Henc, 


4  ~  1/5,  zQ  +  3  «  4a/3,  c  =  a/3,  and 


y  =  31^3(4a)2//3  (z  +  a/3) 1/3 


with  cp'  =  1,  cp"  -  ~  ,  cp" ' 


=  5/da 2  ,  qp'iV  =  -  5 /4a 3  is  indicated 


-....a 


^  ;  ’  .  - .  *>■  .  ^ 


'  'O' 

44-r d,.; ;;;; ;d ;v  <s 


*r  %jt  -s 

XA-XI  - 


li  we  write  tgb  -  j  +  A,  then  we  have  from  Table  5 

and  the  vanishing  of  all  t's  with  3  or  more  subscripts  which 
appear  therein 


D 


~  2  =  -Ay  ,  E 
4a2  “-2 


oa 


oa 


4a 


$ 


2  /I 


'asbb  r  2q  +  ^  +  ( — f~~  +  ~t^)  (~  +  A ) 

ba  4a  a 


2a  -Ay  (Ay  +  —A) 
ooT  8a2  4a2 


~A  A  -  i  A2 
4a  a 


saaabb  *  (~g?)  (|  +  -  2a (1&6 


4" 


64a‘+  32a^  lta 


(§  +  A)  +  4a2  (-I§2  + 

rr  n  ~  O 
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512a  w  256a  I28ac 


&  — i  a2 
a  2a-3  4a2 


5aab  +  3ab  +  sabb 


{a  +  ^  “  |  +  4)  4  4  «  A2 

a 


2"  saabb  +  sab  aanb^ 


1  ,1 


"  a  (a  +  +  +  — Ar)  (~  +  A) 

oa^  24a2  a 


+  3  saaab J 


2a  ( — — ^rr  x  — L_,  J 
64a  *  64a^ 


1  i  1  o 

7  o  -  ~-p  A  -  i-  A41- 
48aJ  3a2  a 


;  r«SH 


T  Y  -  V 3  _ 

-*•***  ~ 


ske  y  =  yh%^  +  7pxB5  *  ycr( 


ji  +  fsjv^*2  +  Tc2)  -  jl  -  ^|uV 


+  Pz  +  a-  +  3 

(4a  aJ  a 


(7A-l)TA3TB2+(rB-a}TB3(TA2+Tc2) 


3/,  2  .  2 


3.  2 


+  (7g-1)Vtb 


J  +  jj  [v^'aV^c^'cV 


r  ii2  ■  ^I(^3tb+vb3^a%2)^c3^) 

("  7t2 + aI^v+vaV) 


c 

+  terms  of  order  >  t 


elo  y  -  Z(ra-3)xa4 


+  terms  of  order  >  t^, 


If  we  put  a  -  10,  then  the  variance  formula  becomes, 
to  five  decimals, 


1X^1 


4^ 


+  c.oooio  ig2 ( 

-  0.00002  (7ATAV+>BTB3(TA2+Tc2)+r0rc3TA2) 

-  0.00002  (yA ta3t02  +  rcTA2Tc3) 


+  .  .  . 

and  clearly  the  correction  terms  are  all  very  negligible. 
Before  transformation,  we  had 

*<■ 

tAB  =  fcBC  =  if?  =  S  +  h  =  °-11- 
<=p  a 

fcAC  =  f?  =  I  -  3  =  °-°99  , 

P 

and  the  variance  formula  was 

var  z  -  V  +  tb2  +  Tc2 

+  °-01210  ^B2  W2  +  V2 

+  0,00980  rA2Tc2 

Thus  we  have  reduced  the  first  correction  term  by  a  factor  of 
100  or  so,  although  it  was  already  rather  small. 


.  . 


Tilts-  above  rtsuiy  were  cased  on  the  plausible* 


seeming  assumption  that  the  leading  terns  in  qp"  and  qp'M  suffice 
to  select  an  adequate  transformation.  It  will  prove  enlighten¬ 
ing  to  Inquire  into  what  the  "best"  transformation  proves  to  be  for 
a  =  10.  We  find,  using  all  terms,  <p"  =  -.05225,^"'  =  .004  870  03125 


1  -  p 


.002  730  062  5 
7602  139  968 


1.2 757+ 


zQ  +  c  =  24.416 


24.  ^e1*2^ 

- . 275749 


-213.705 


j 


y  =  -213.705  (z+l4.4l6)*-2757 


Tha.s  transformation  Is  clearly  quite  different  In  appearance 
from  the  previous  one,  although  Its  derivative  values  are  quite 
similar,  namely 


ii 

Approximate 

Exact 

T 

- . 0500 

-.05225 

All 

4> 

. 00625 

. 00487 

91V 

-.00125 

-.OOO65 

9V 

. 00034 

.00011 

X  GLOSSARIES  AND  NOTATION 


i 


* 


I 


' 1 


36.  Pie  X*  notation. 

We  have  used  at  various  points  the  very  convenient,  hut 
not  standard,  Z*  notation.  With  its  convenience  there  comes  a 
need  for  some  care  in  its  use.  At  this  point,  after  recalling  its 
definition,  some  examples  will  be  provided.  In  general,  2*  implies 
the  summation  of  all  the  different  terms  of  the  form  written  after 
it  which  can  be  obtained  without  identifying  subscripts.  In  this, 
its  usage  corresponds  to  our  usual  practice  in  writing  down  simple 
formulas.  We  write 

p  p  p 

(x1+x2)  =  x]_  4:'2  +2X1X2 


but 


(x1+x2)(y:i+y2)  =  x1y1+x2y2+x1y2+x2yi 

Similarly  we  now  write 

( 2x  ) 2  =  Lx,  2+22*xx, 

'a '  a  a  b 


but 

(Sxa)(Syb)  =  Zxaya+Z*xayb . 

#  # 

Note  the  22  in  one  formula  and  the  Z  without  the  two  In  the 

p 

other.  This  exactly  corresponds  to  our  writing  2x^x2  in  (x^+x^)- 
and  l(x1y2+x2y1)  in  (x1  x2)(y1+  y2) . 

As  a  consequence  we  must  be  careful  of  what  happens 

when  we  specialize  variables,  break  up  terms,  or  multiply  expres- 

*  * 

sions.  Thus  if  we  put  y^x^  ln  X  x&yb  the  answer  is  22  x&x^  and 


~  ' 


not  merely  Z  As  another  example,  conr 1 der  Che  following 


*  *  t  . 

Z  x  y.  =  I  ( x  v, -x  y  +x,v  ) 

ci  k  a.  O  a.  at  d  a 


“ 2  (xa<v  ya> +  v.) 


“  1  xa(yb-ya)  +  (n-1)  aaya 


where  we  find  n-1  times  Zx  v  =  Z~x,,y„  rather  than  merely  Z  x  y  . 

a  a  a  a  a  a 

(Here  n  =  the  number  of  values  taken  on  by  a.)  Finally,  consider 

.  #  p  p 

(Z  x  xfe)  ,  which  we  now  write 


(Z  xa2xb)(Z  xc2x<j)  =  2*xa^xb2 


*-»  xv  -A,-. 

a  b 

#-3  0 

+2  Z  x0Jx,  J 
a  b 

#  4 

+2  Z  x_  x.x , 
a  b  d 

#  2  ^ 

+2  Z  xa  xb^xd 

..#222 
+6  Z  xa  xfc  xc 

*  2  2 
+  Z  xa  xfexc  xd  . 


It  would  not  be  too  difficult  to  miss  some  of  these  numerical 
coefficients . 

All  these  needs  for  care  considered,  however,  the  ad- 

# 

vantages  of  the  Z  notation  seem  to  outweigh  the  disadvantages. 
39 .  Glossary  of  statistical  and  other  special  terms  used . 
average  Arithmetic  mean,  especially  of  a  probability  distribu¬ 
tion.  (within  range  of  standard  usage.) 
average  point  Hie  situation  in  which  each  individual  variable 

takes  on  its  average  value.  (special) 


** **‘***’*<Si^^ 


X-j 


Cocunulant  A  sanlnvari  ant-  ~ .  .. 

-  ~ant  measure  of  mutual  variation.  Defini¬ 

tions  of  ,he  first  thre-*  and  statement  of  their 

properties  given  In  Section  19.  (new,  not  yet 
standard) 

i^£JS2S2tlor  The  fourth  cocumulant  -  a  function  of  four  argu¬ 
ments.  (See  Section  19  for  definition  and  proper- 
(new,  not  yet  standard) 

correlated  Not  having  zero  covariance,  (standard) 
fewness  The  third  coolant  -  a  function  of  three  argu¬ 
ments.  See  Section  19  for  definition  and  proper- 
ties.  (new,  not  yet  standard) 

covariance  The  average  proiuct  of  deviations  from  means  -  the 

second  eocumulant .  f  a  Section  19  for  definition 
and  properties,  (standard) 

cumuli ^.^ejnl^rlantl  A  seminvarlant  polynomial  In  the  mo¬ 
ments  of  a  distribution;  a  coefficient 
in  the  expansion  of  the  logarithm  of 
the  moment  generating  function  in  terms 

-i  k 

01  f  A ••  The  first  cumulants  are 
average,  variance,  skewness  and  elon¬ 
gation.  (standard) 

Function  expressing  a  normal  die- 
trlbution  In  terms  of  the  total 
probability  of  values  ]ess  than 
any  given  value,  (standard) 


— A  quantity  with  a  probabillcy  distributer  hr-r  » 

(standard)  *  16  “«• 

— p ...r  1 ra ^  '  v e  vaJ  U  j  A  valup  or  a  dpr-fvn*-*v  ♦ 

- -  Ot  a  a.riva^w  at  a  preassigned  point. 

(not  yet  widely  used) 

elongation  The  fourth  cumulant  usual! v  «  m^or. 

'  USjally  a  measure  of  longtailedness . 

See  also  section  4.  (new) 

independence  (statistical)  Two  chance 

- - 1  quantities  are  statistically  independent 

Whan  Knowledge  of  the  value  of  one  gives  no  Information 

about  the  probability  distribution  of  the  other. 

( standard) 

—  ^  ™1Ue  —  P0«r  of  the  quantity  conceded 

See  also  Sections  a  to  It.  (standard) 

2°5S5L£i£trlbution  A  particularly  simple  shape  of  distribution  of  probability 

which  can  be  characterized  in  mny  ways  ;  as  the  distribution 
of  the  sum  of  an  Indefinitely  large  number  of  independent 
quantities,  as  a  distribution  all  of  whose  higher  cumulants 
vanish,  as  an  example  of  the  normal  "bell -shaped  curve". 

AlBO  "Gaussian"  or  "Maxwelllon" .  (standard) 

ES__a_effect  The  effect  of  varying  one  individual  variable,  while  holding 
all  the  others  constant,  (not  yet  standard) 
relative  derivative  values  Katies  of  higher  derivative  values  to  the 

corresponding  products  of  first  derivative  values . 

(not  yet  standard) 


-  x-5  - 


relative  momenta  Ratios  of  higher  moments  about  the  mean  to  that 

power  of  the  standard  deviation  which  makes  the 
result  dimensionless.  (not  infrequent) 
response  A  variable  whose  value  is  regarded  as  caused  by  the 

values  of  another  variable  (stimulus)  or  other  variables 
(Individual  variables,  stimuli).  (standard) 
semlnvarlant  (see  cumulant) 

skewness  The  third  cumulant  (=  the  third  of  moment  about  the  mean), 
usually  a  measure  of  asymmetry  In  long-taiiedness .  See 
also  Section  4.  (marginal) 

standard  measure  A  distribution  is  expressed  In  standard  measure 

when  its  average  is  zero  and  its  variance  unity. 

( standard) 

standardized  deviate  A  quantity  with  average  zero  and  unit  vari¬ 
ance.  (standard) 

terms  The  type  of  unit  in  which  a  variable  is  expressed  to  be  con¬ 
trasted  with  "scale"  which  Indicates  the  size  of  the  unit 
used.  (The  variable  "temperature"  for  example,  can  be  ex¬ 
pressed  In  terms  of  degrees,  log  degrees,  reciprocal  de¬ 
grees,  etc.)  (new-~no  clear  pattern  of  usage) 
variance  The  root-mean-square  deviation  from  the  average;  the 
second  cumulant.  See  also  Section  2.  (standard) 

40 .  Glossary  of  abbreviations. 

"ave"  Average  value  of  expression  following.  (Used  by  some 
authors . ) 

"coe(^, -,-,-)"  coelongation  of  four  expressions  appearing  as 

arguments  (not  yet  standard) 


" cok ~ ) "  coskewn&ss  of  three  ex&r essions  appearing  as  arguments. 

(not  yet  standard) 

Hcov(-,-)"  covariance  of  two  expressions  appearing  as  arguments 
( standard) 

"elo"  elongation  (fourth  cumulant)  of  expression  following. 

(not  yet  standard) 

"ske"  skewness  (third  cumulanr)  of  expression  following,  (not 
yet  standard) 

"var"  variance  (second  cumulant)  of  expression  following. 

( standard) 

4l.  Notation  used  here  for  response  functions. 

The  main  topic  of  this  memorandum  deals  with  the  behavior 
of  a  response  as  a  function  of  individual  variables.  The  following 
notations  are  used  consistently 

(a)  No  special  assumptions: 

z  -  h(w1,w2.,  . .  .  ,wk) 

(b)  Higher  unmixed  derivatives  all  vanish: 

z  =  g(vi> v2>  *  *  ”vk) 


(c)  Transformed  response,  with  individual  variables  such  thac 
all  higher  unmixed  derivatives  all  vanish: 

y  =  f (ux,u2, . . . ,uk) 


in  addition,  various  other  notations  are  use 
eluding 


id  transiently,  in- 


< 


y 


9 


2 


z  =  h(wa,wb,wc), 
z  =  h(wA,wB,wc), 
z  =  f(VA.VB#Vc), 


with  restrictions  stated  in  context. 

42.  Index  of  notations  used  "on  the  line" 

The  index  to  symbols  which  follows  is  divided  into 
three  sections:  first,  in  this  section,  symbols  usod  "on  the 
line",  next,  in  the  next  section,  symobls  used  as  subscripts, 
and  finally,  in  section  44,  symbols  used  as  exponents  on  other 
superscripts.  Uses  of  limited  extent  are  specified  as  inciden¬ 
tal  . 

Usage  "on  the  line" 

al>a2, . . . ,ak  Incidental  constants  (abstract  only.) 

A  A  choosable  constant  (sections  16,  33,  37) 
b  An  incidental  constant  (section  9  only) 

c  An  incidental  constant  (section  9  only),  a  choosable  constant 
(section  16,  33>  37) 

C'°  V°21*°a*  etc.  Capacitance  of  element  indicated  by  subscript 

(in  delay  line  example),  (section  lOff) 


D  Abbreviation  for  <pf *  '  W.  (section  1 6  only) 

K  Abbreviation  Tor  sp"1?"*  (section  15  only) 

f  ^  function  of  k  arguments.  (sections  4l,  29,  15ff) 

^a' ^bc' ^3* ^112'  e“c  •  Values  of  the  Indicated  derivatives  of 

f  (-,-,...,-)  at  the  average  point  --  i.e., 
with  each  argument  of  f (-,-,...,-)  at  its 
average  value.  (sections  15ff) 
s( -> ->•*., -)  A  function  of  k  arguments,  (sections  4l,  29, 

various) 


sa,sbc's3'sli2*  etc*  Derivative  values  of  g(~, at  the 

average  point  (cp.f&,  etc.  above) ( sections 
7ff) 

^ (-»-»•••» -)  A  specific  incidental  function  of  k  arguments. 

(section  7  only) 

G, G&  Relative  fifth  moments  about  the  average,  (sections  2,3,5 
and  21-23; . 

-j . . . ^ -)  A  function  of  k  arguments,  (section  4l)  29,  various) 
ha,hbc-'h3,hll2'  etc*  Derivative  values  of  h(-, at  the 

average  point  (cp.f  etc.  above), 

(sections  3ff) 

H( -,-,-)  Attenuation  as  function  of  WA,Wfi, Wc  (section  27  only) 
j  Number  of  sections  of  delay  line  (sections  10,  11,  24ff ) . 
k  Number  of  arguments  in  response  functions  =  number  of  indi¬ 
vidual  variables  treated  as  affecting  response,  (sections 
4l,  various) 


L>L'..^2i- l'La>  etc  *  Inductance  of  element  indicated  by  subscript 

(in  delay  line  example .  (section*  10,  11, 

24ff .) 

p  See  exponent  In  section  44. 

p(z)dz  Generic  probability  density  of  z.  (section  2  only) 
qa  Incidental  deviation  of  w&  from  its  average,  (section  23 
only) 

qat),  etc.  Incidental  expressions  (section  33  only) 

Ra,RB'RC  Values  resistance  in  attenuator  example  (sections 
2,  27ff) 

sab,saab*  etc*  Relative  derivative  values  of  f,  e.g.  a  ^  = 

fat/fafb'  (sections  15ff,  33ff) 

tab'^aab'  etc*  Relatlve  derivative  values  of  g,  e.g.  t  = 

gftv/g-gh*  (section  8ff) 
a  d  u 

ul,u2'***'uk  A  set  of  individual  (=  component)  variables). 

(sections  4l,  29,  15ff) 

ua,ub'***  General  examples  of  u1,u2,„..,uk.  (various) 
v  A  chance  quantity  (Part  V  only) 

v1,Vg, ...,vk  A  set  of  individual  (=  component)  variables  (sections 

41,  29,  6ff) 

vl^'v2^,  .  •  •^vk^1^  An  incidental  set  of  particular  individual 

(=  component)  variables  (section  7  only) 
va(ua)  A  change  in  terms  of  the  ath  individual  variablj  (sections 
15,  29,  30) 


""•r*  ,u 


-  x-xo  - 


vi  }X2  *  "  *,vk  ^hddent&l  notation  for  the  average  values  of 

V1  * v  2*  ‘  vk'  (se^tion  9  only) 

Incidental  notation  for  the  coordinates  of  the 

XX  Ki 

new  average  point  (section  9  only). 

v2i-J,v21  Measures  of n/^21-1  and \/C2i  respectively  (sections 
25,  26  with  differing  usage.) 

v  I  il  v  v  iv,  etc.  Successive  derivative  values  of  v  (u  ) 
aa'a'a'  a'  a' 

at  the  average  point.  (sections  30ff.) 
VA7VB,VC  Measures  of  deviation  of  input  shunt  conductance,  series 
resistance  and  output  shunt  conductance,  respectively, 
from  their  nominal  values  (in  the  attenuator  example) 
(sections  12,  28) 

VA'VB’VC  Other  similar  measures  (section  28  only) 

w-,,w0, - ,vr,  A  set  of  k  Individual  (=  component)  variables 

i  2  k 

(sections  4l,  29,  various) 

w.,w-,w~  Measures  of  deviation  of  element  resistances  from  nominal 
A7  B7  C 

values  (in  attenuator  example) ( sections  12,  27) 

WA,VfB;,Wc  Similar  measures  (section  27  only) 

x  A  chance  quantity  (Part  V  only) 

Xi,x2,...,x  A  set  of  k  individual  (=  component)  variables. 

(sections  41,  various.) 

y  A  response  variable  (sections  4l,  various)  A  chance  quantity 

s. 

(Part  V  only) . 

Ya,Yc  Conductances  of  shunt  elements  (in  attenuator  example), 
(sections  12,  28) . 


*  X«  J.1  « 


m 


i 


r 


2  A  response  variable, 

(Fart  V  only) 


{ sections  *♦!,  various) 


A  chtne*  quantity. 


a  Nominal  attenuation  fin  thp  v  ✓ 

in  me  attenuator  example)  (Sections  12, 

27,  36,  37) 

@1'P2  Dimensionless  functions  of  higher  moments .  (Standard. 

Sections  4,  14) 

P,r,6  Specific  rational  functions  of  a  (in  the  attenuator  example* 
(sections  12,  27,  28,  36). 

W  etc-  Relative  third  moments  (in  the  attenuator  example 
tnls  meaning  applies  only  to  7-3  with  subscripts.) 
(Moderately  used,  section  2ff.) 

7V72  Dimensionless  functions  of  higher  cumulants  (standard. 
Sections  4,  13.) 

T,ra  Relative  fourth  moments.  (Transient,  Section  2ff.) 

A  Deviation  of  t&b  from  1/a  (incidental),  (section  26  only) 

S,5*  "small"  deviations  from  nominal  specifying  upper  tolerances 
(sections  l,  3,  4) 

e,e*  "small"  deviations  from  nominal  specifying  lower  tolerances 
(section  1) 

^21-1*^21  Coefficients  of  variation  of  and 

respectively  (in  delay  line  example).  (sections  11, 

26) 

Ki  lth  cumulant  =  lth  seminvariant.  (standard,  section  4  only) 

^1  lth  moment  about  the  average  [about  the  mean],  (standard, 
sections  2,  4) 

t 

M-i  lth  moment  about  the  origin  (standard,  section  2  only) 


Z~t2  - 


vl,v2*  *  *  vk  W*p2*c««eBt  of 

*  -  *5®  ( *n  response  imitc), 

(section  9  only) 

<?(-)  Function  expressing  a  transformation  of  response  (section 

25,  16,  2jTT . ) 

<r.  *■  i  V 

'9  *9  ,9  ,  etc.  Successive  derivative  values  of  <p(.)  at  the 

new  average  point  (sections  15,  16,  3Qff.) 

^'a(wa)  Incidental  change  In  terms  of  the  ath  Individual  variable, 
(section  7  only) 

2*  sign  of  restricted  summation.  (Discussed  In  section  38) 
o,cra,  etc.  Standard  deviation  =  square  root  of  variance  =  root- 
mean-square  deviation  (standard,  section  2ff) 

2  Sign  of  unrestricted  summation,  (standard,  various) 

Ta  RMS  value  of  Individual  contribution  to  response  -  t  a  • 
fa(var  xa)1/2-  (Section  8ff.)  3  3 

JUu — Index  of  notations  used  as  subscripts . 

a  Generic  subscript,  usually  ranging  from  1  to  k. 

A  Identifies  input  shunt  element  In  attenuator  example 
(sections  12,  27  ff) 

b  Llke  a'  but  ordinarily  distinct,  (various) 

B  Identifies  series  element  In  attenuator  example 
(sections  12,  27 ff) 

0  Like  a,  but  ordinarily  distinct  from  both  a  and  b. 

(various) 

°'  IdenUfleS  0utput  shunt  el®ent  in  the  attenuator  example, 
(sections  12,  27ff) 

1  Index  ranging  from  1  to  j.  (sections  10,  11,  24ff) 


-  j t-ii  - 

0  Corresponding  to  average  point,  (various) 

o  Corresponding  to  average  point .  (abstract 

44.  Index  of  notation  used  as  superscripts. 

i,j  Incidental  exponents,  (section  21  only) 

m,n  Incidental  exponents,  (section  21  only) 

p  Exponent  in  power  transformation,  (sections  16,  33.  37) 

p,q  Incidental  exponents,  (section  21  only) 


